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developed  a  new  method  of  synthesis  via  the  nitrodephosphorylation  of 
phosphoramides ,  which  is  very  convenient  for  certain  nitramines. 
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INTRODUCTION 


The  principal  objective  of  this  project  has  been  to  identify  and  deter¬ 
mine  the  thermal  rate  parameters  for  the  initial  and/or  rate-limiting  steps  in 
the  gas -phase  decomposition  of  nitramines.  In  support  of  this  goal,  which  has 
also  received  funding  from  the  Army  Research  Office  (ARO) ,  we  have  applied 
three  different  approaches.  Experimentally,  we  have  used  two  laser  pyrolysis 
systems  to  study  the  initial  and  secondary  reactions  in  the  decomposition  of 
dime thy Ini tramine .  The  theoretical  effort  has  involved  the  use  of  ab  initio 
quantum  mechanical  calculations  to  help  determine  the  extent  to  which  nitro- 
nitrite  rearrangement  is  a  competitive  pathway  in  the  decomposition  of  nitra¬ 
mines.  The  results  of  each  of  these  approaches  are  summarized  in  this  section 
and  presented  in  detail  in  Appendices  A  through  D,  each  of  which  is  a 
manuscript  submitted,  or  to  be  submitted,  for  publication. 

As  an  important  class  of  energetic  materials,  nitramines  have  been  the 
subject  of  numerous  studies  whose  goals  have  been  to  elucidate  various  facets 
of  the  chemistry  that  determines  whether  a  self-sustaining  reaction  is 
achieved  and,  in  the  case  of  propellant  compositions,  determines  the  rate  of 
combustion.  Inasmuch  as  attempts  to  modify  the  burn  rates  of  nitramine 
propellants  have  met  with  only  very  modest  success,  and  inasmuch  as  the  high 
pressure -dependence  of  HMX  combustion  rate  aggravates  a  tendency  towards 
instability  and  deflagration-to-detonation  phenomena,  achieving  a  better 
understanding  of  nitramine  chemistry  could  have  substantial  practical 
impact.  However,  the  chemical  details  of  nitramine  decomposition,  like  those 
for  other  energetic  materials,  have  been  very  difficult  to  ascertain.  This  is 
hardly  surprising,  because  in  substances  that  decompose  by  one  or  more  rate- 
limiting  initial  endothermic  steps  followed  by  a  sequence  of  highly  exothermic 
and  very  rapid  steps,  the  nature  of  the  rate-limiting  steps  themselves  tends 
to  be  obscured.  Characterizing  nitramine  decomposition  is  particularly 
difficult  not  only  because  secondary  bimolecular  reactions  of  radicals  with 
NC>2  are  particularly  fast  but  also  because  the  aliphatic  radicals  formed  from 


the  carbon-nitrogen  skeleton  have  a  w^der  range  of  possible  reactions  than, 
for  example,  the  typical  aromatic  radical  produced  in  nitroaromatic 
decomposition. 


Dimethylnitramine  (DMNA)  decomposition  is  not  exactly  analogous  to 

decomposition  of  the  cyclic  nitramines,  which  have  several  additional  pathways 

1  2 

potentially  available.  ’  However,  because  previous  research  shows  that 

1  Q 

understanding  DMNA  decomposition  is  difficult,  we  chose  to  study 
competition  in  the  reaction  pathways  of  DMNA  before  tackling  RDX  and  HMX, 
which  have  more  complex  chemistry,  lower  vapor  pressures,  and  higher 
reactivity. 

The  expedient  commonly  used  in  attempts  to  minimize  the  obscuring  of 
rate-limiting  steps  by  rapid  secondary  reactions  is  to  perform  the  decomposi¬ 
tion  at  much  lower  temperatures  and/or  much  lower  pressures  (concentrations) 
than  those  involved  during  "normal"  operation.  Drastically  different 
temperature  and  pressure  conditions  are  likely  to  cause  significant  shifts  in 
competing  reactions  and,  therefore,  leave  undetermined  the  rate-limiting  steps 
that  are  actually  operative  under  practical  conditions.  In  our  research,  we 
lowered  the  pressure  but  not  the  temperature  and  used  very  short  reaction 
times . 


LASER  PYROLYSIS  OF  DIMETHYLNITRAMINE  AND  DIMETHYLNITROSAMINE 


To  help  determine  the  rate- limiting  step  or  steps  at  relevant 
temperatures,  we  chose  an  intermediate  experimental  course  in  which  the 
pressure  but  not  the  temperature  is  significantly  lowered  and  the  reaction 
time  is  very  short.  A  pulsed  infrared  laser  is  used  to  heat  the  substrate  by 
heating  an  absorbing  but  unreactive  gas  (e.g.,  SFg)  that  then  transfers  its 
kinetic  energy  by  collision;  in  this  way,  the  substrate  is  brought  to  high 
temperature  (700-1200K)  in  1-3  fxs .  Reaction  then  occurs  in  the  1  to  15  fis  in 
which  the  acoustic  expansion  wave  moves  radially  inward  toward  the  center  of  the 
laser-heated  region.  By  using  different  methods  of  sampling  the  product 
mixture,  we  can  examine  different  extents  of  secondary  reactions. 

We  have  used  two  different  techniques  to  determine  products:  One  involves 
passing  the  contents  of  the  cell  shown  in  Figure  1  through  an  automatically 
operated  gas -sampling  valve,  which  periodically  injects  a  sample  into  a  GC/MS 
system.  This  provides  a  time-averaged  mixture  of  products  from  the  entire 
reaction  zone  and  a  reliable  average  temperature/reaction  time  from  the 
simultaneously  determined  fractional  decomposition  of  a  temperature 
standard. The  second  technique  involves  a  separate  apparatus  (Figure  2) 
in  which  the  laser-heated  region  is  located  immediately  adjacent  to  the 
sampling  nozzle  of  a  molecular-beam  mass  spectrometric  detection  system.  In 
this  case,  the  gas  mixture  in  the  laser-heated  region  expands  toward  the 
nozzle,  and  the  product  mixture  is  frozen  by  the  molecular-beam  formation.  The 
product  mixture  in  the  periphery  of  the  laser-heated  region  has  the  least 
opportunity  for  secondary  reaction,  and  that  in  the  inner  portions  has 
progessively  more  opportunity.  With  both  sampling  techniques,  the  extent  of 
secondary  reaction  can  be  varied  by  varying  the  substrate  partial  pressure,  the 
extent  of  the  initial  reaction,  and  the  presence  of  radical  scavengers  and 
radical  traps.  The  presence  of  scavengers  and  traps  has  been  particularly 
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GC/MS -Monitf'-ed  Laser  Pyrolysis 


DMNA  was  pyrolyzed  over  the  temperature  range  800  to  900K  in  a  flow  system 
using  pulsed  infrared  laser  heating.  The  major  product  under  all  conditions 
was  dimethylnitrosamine  (DMNO) ,  with  smaller  amounts  of  N-methylmethyleneimine 
and  dimethylamine  also  produced  when  radical  scavengers  and  traps  were  not  both 
present.  Arrhenius  parameters  for  the  decomposition  were  log  k(s’^)  -  13.5  + 
0.6  -  (37.4  +  2.5)/2.3RT.  These  parameters,  which  were  produced  under  a  wide 
range  of  conditions  (different  scavengers,  different  temperatures,  and  varying 
amounts  of  added  NO  as  a  trap  for  dimethylamino  radicals),  are  not  consistent 
with  simple  bond  scission  as  the  sole  rate- limiting  step.  The  principal 
product,  DMNO,  was  studied  under  the  same  conditions;  DMNO  has  no  plausible 
decomposition  pathways  other  than  simple  bond  scission,  and  it  decomposed  with 
parameters  appropriate  to  such  a  process,  namely  log  k(s'^)  -  15.8  ±  1.1  - 
(50.3  +  3.4)/2.3RT.  Moreover,  the  temperature  dependence  was  evidently  not 
distorted  (as  feared)  by  secondary  reactions  involving  the  highly  reactive 
methyleneimine ,  which  was  the  major  product  of  DMNO  decomposition,  but  a  less 
important  product  of  DMNA  decomposition.  The  experimental  products  and 
temperature  dependence  can  be  reproduced  via  a  mechanistic  numerical  model  when 
N-NO2  bond  scission  and  nitro-nitrite  rearrangement  are  competitive  initial 
steps  and  the  displacement  of  NO2  from  DMNA  by  NO  is  included  as  a  low- 
temperature  route  to  DMNO.  Thus,  these  results  lead  us  to  conclude  that  DMNA 
decomposes  not  only  by  the  expected  simple  bond  scission  reaction  but  also  by 
an  intramolecular  nitro-nitrite  rearrangement.  Such  reactions  have  been  shown 
to  occur  in  C-NO2  systems  but  have  not  been  previously  reported  in 

nitramines.  These  results  are  described  in  detail  in  Appendix  A,  a  manuscript 
submitted  to  the  Journal  of  Physical  Chemistry. 


Molecular -Beam  Mass-Spectrometrically-Monitored  Laser  Pyrolysis 

Because  the  products  of  an  initial  nitro-nitrite  rearrangement  can  also 
result  from  a  secondary  bimolecular  reaction  in  which  the  dimethylamino  radical 
is  oxidized  by  NO2,  and  because  the  evidence  summarized  above  for  the 
intramolecular  rearrangement,  while  extensive,  is  largely  indirect,  we  used  a 
second  laser-pyrolysis  system  to  examine  directly  the  products  of 
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dime thy Ini tramine  decomposition.  With  this  system  (Figure  2) ,  the  product 
mixture  in  the  laser-heated  region  expands  directly  into  the  nozzle  of  a 
molecular-beam  sampling  system  and  is  detected  by  a  quadrupole  mass 
spectrometer.  The  initial  products  include  not  only  dimethylamino  radical  and 
NO2  but  also  comparable  amounts  of  NO  and  various  fragment  ions  of  the  nitroxyl 
radical  (CH^^NNO-.  Although  the  bimolecular  oxidation  of  radicals  by  NO2  is 
known  to  be  very  rapid,  the  NO2  partial  pressure  is  too  low  and  the  time 
available  before  the  product  mixture  is  frozen  in  the  expansion  too  short  for 
this  route  to  be  important.  In  fact,  when  the  partial  pressures  of  initially- 
produced  NO2  and  dimethylamino  radicals  are  increased  by  about  a  factor  of  ten, 
the  consumption  of  NO2  in  secondary  bimolecular  reactions  can  be  observed. 
Calibration  of  the  detection  system  using  known  flows  of  laser-heated  NO  and 
NO2  indicates  the  rearrangement/N-N02  bond  scission  branching  ratio  (at  about 
900K)  is  -  0.7.  Thus,  these  results  confirm  the  inferences  drawn  from  the  low 
temperature  dependence  consistently  observed  in  the  GC/MS -monitored  nitramine 
pyrolysis.  This  work  is  described  in  detail  in  Appendix  B,  a  manuscript  to  be 
submitted  to  the  Journal  of  Physical  Chemistry. 

In  ongoing  studies  under  ARO  support,  we  are  probing  the  laser  pyrolysis 
reactor  still  more  directly,  using  laser- induced  fluorescence  (LIF)  to 
determine  the  nascent  N0/N00  ratio. 


AB  INITIO  QUANTUM  MECHANICAL  CALCULATIONS 


To  further  establish  the  significance  of  nitro-nitrite  rearrangement  as  an 
initial  step  in  decomposition  of  DMNA,  for  which  we  have  obtained  evidence  in 
both  of  the  laser-pyrolysis  systems,  we  have  undertaken  a  theoretical  study  of 
the  rearrangement  process  in  the  analogous  unsubstituted  nitramine.  Structures 
of  NH2NO2 ,  the  four  isomers  of  NH2ONO,  dissociation  products,  and  the 
transition  state  have  been  optimized  using  MCSCF  analytic  second  derivatives 
with  a  4-31G  basis  set.  The  geometry  of  the  transition  state  may  be  described 
as  separated  NH2  and  NO2  fragments  held  together  by  dipole-dipole  interaction. 
The  N-N  bond  length  is  2.85  A,  twice  as  long  as  at  the  NH2NO2  equilibrium 
geometry.  The  wave  function  of  the  transition  state  is  characterized  by  two 
singly-occupied  orbitals.  Optimization  of  the  transition  state  depends, 
therefore,  on  the  use  of  MCSCF-based  methods. 

Energies  have  been  obtained  by  large-scale  multireference  single-  and 
double-excitation  Cl  calculations  with  a  6-31G  basis  set.  The  transition 
state  is  predicted  to  lie  46.68  kcal/m  above  OT^NOj  and  0.87  kcal/m  below  the 
NH2  +  NO2  asymptote.  Because  the  significant  frequencies  of  the  transition 
state  are  very  similar  to  those  of  the  separated  fragments,  the  A  factor  for 
rearrangment  should  be  not  much  less  than  that  for  bond  scission.  Assuming 
that  results  for  the  unsubstituted  species  may  be  applied  directly  to  dimethyl - 
nitramine,  the  rates  for  nitro-nitrite  rearrangement  and  N-NO2  bond  scission 
should  be  quite  comparable  in  the  latter  system,  in  agreement  with  the  results 
from  the  molecular-beam  mass-spectrometrically  monitored  laser  pyrolysis 
experiments  reported  above.  However,  we  note  that  quantitative  agreement  with 
the  temperature  dependence  observed  in  the  GC/MS-monitored  LPHP  system  would 
require  an  A- factor  and  activation  energy  substantially  lower  for  the 
isomerization  in  order  that  the  overall  temperature  dependence  would  be  roughly 
midway  between  that  for  each  of  the  competing  initial  steps. 

Details  of  the  theoretical  study,  which  has  been  performed  in  collabora¬ 
tion  with  Dr.  Megumu  Yoshimine  of  IBM  Almaden  Research  Center,  are  presented  in 
Appendix  C,  a  draft  of  a  manuscript  to  be  submitted  to  the  Journal  of  Physical 
Chemistry. 


NITRAMINE  SYNTHESES 


In  the  course  of  repeated  syntheses  of  dimethylnitramine,  we  developed  a 
new  route  that  is  extremely  convenient  and  safe  for  synthesizing 
dimethylnitramine  and  presumably  other  nitramines  as  well.  This  method 
involves  nitration  of  the  respective  phosphoric  triamide  and  avoids  the 
nitrations  that  can  be  very  hazardous,  even  at  moderate  scale.  For 
dimethylnitramine,  this  method  is  extremely  cheap  and  convenient,  because  it 
simply  requires  the  commonly  used  solvent  hexamethylphosphoramide .  The 
synthesis  of  dimethylnitramine  from  hexamethylphosphoramide  is  described  in  the 
short  manuscript  included  as  Appendix  D. 
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CONCLUSIONS 


The  result  of  the  three  approaches  to  understanding  dimethylnitramine 
decomposition  we  have  used  in  this  research  is  a  demonstration  of  the  existence 
of  a  new  initial  reaction  pathway  for  nitramines,  a  unimolecular  nitro-nitrite 
rearrangement,  followed  by  rapid  scission  of  the  extremely  weak  (Me^NO-NO 
bond.  This  pathway,  which  tends  to  be  competitive  with  the  expected  N-NO2 
scission,  is  analogous  to  that  shown  recently  for  aliphatic-  and  aromatic- 
nitrocompounds,  but  had  not  been  previously  reported  for  nitramines.  Recent 
dimethylnitramine  decomposition  results  in  the  literature  had  in  fact,  been 
interpreted  as  evidence  for  exclusive  N-NO2  fission.  The  new  nitro-nitrite 
rearrangement  pathway  may  have  relevance  for  the  practical  functioning  of 
nitramines,  since  it  leads  directly  to  a  more  oxidized  skeleton  (and  possibly 
to  the  various  RDX  nitrosamine  products  recently  reported  by  Hofsommer  and 
Glover  in  condensed  phase  decomposition^),  but  may  well  block  the  concerted 
(or  sequential)  depolymerization  to  the  nitro-imine  monomer  units.  The 
importance  of  the  latter  products  is  that  they  can  lead  directly  to  ^0  and 
formaldehyde,  and  thence,  in  a  highly  exothermic  reaction,  to  N2,  CO,  and 
H20.  Thus,  changes  in  the  branching  ratios  between  nitro  nitrit »  rearrangement 
and  competing  initial  steps  (N-NO2  scissions  and,  in  cyclic  nitramines, 
depolymerization)  may  well  help  determine  whether  the  decomposition  process 
proceeds  to  a  self  sustaining  stage.  In  more  general  terms,  the  results  serve 
to  emphasize  the  multiplicity  of  pathways  available  to  nitramines,  and  the 
necessity,  difficulties  notwithstanding,  for  characterizing  these  pathways  if 
the  factors  controlling  nitramine  ignition,  initiation,  and  deflagration- to - 
detonation  behavior  are  to  be  really  understood. 
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ABSTRACT 

Pyrolysis  of  dimethylnitramine  (DMNA)  and  dimethylnitrosamine  (DMNO)  was 
carried  out  in  a  flow  system  over  the  temperature  range  800  to  900  K  using 
pulsed  infrared  laser  heating,  via  SFg  in  a  250  torr  bath  of  COj  and  radical 
scavenger.  Arrhenius  parameters  for  DMNO  and  DMNA  decomposition  were 
log  k(s_1)  -  15.8  ±  1.1  -  (50.0  ±  3.4)/2.3RT,  and  13.5  ±  0.6  - 
(37.4  ±  2.5)/2.3RT,  respectively.  The  former  set  of  parameters  is  consistent 
with  simple  bond  scission  as  the  rate  limiting  step;  the  latter  set,  which  was 
produced  with  different  scavengers,  temperature  standards,  and  varying  amounts 
of  added  NO  as  a  radical  trap,  is  not  consistent  with  simple  bond  scission. 

The  experimental  results  can  be  reproduced  via  a  mechanistic  numerical  model 
in  which  N-NO2  bond  scission  and  nitro -nitrite  rearrangement  are  competitive 
initial  steps  and  the  displacement  of  N0£  from  DMNA  by  NO  is  included  as  a  low 
temperature  route  to  DMNO. 
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INTRODUCTION 
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Some  years  ago  Fluornoy1  reported  that  the  decomposition  of 

N,N-dimethylnitramine  in  static  bulb  experiments  over  the  temperature  range 

165  to  200  K  resulted  in  formation  of  N,N-dimethylnitrosamine  in  roughly  85% 

yield.  He  explained  this  product  via  recombination  of  NO  and  dimethylamino 

radicals.  This  reaction  is  in  itself  completely  reasonable.  However, 

Fluornoy  found  it  necessary  to  invoke  reactions  that  would  produce  NO  in  high 

o 

yield,  but  these  are  somewhat  vague  and  questionable,  as  discussed  by  Lin  . 
Several  years  ago,  weJ  used  the  very  low-pressure  pyrolysis  technique  (VLPP) 
to  examine  dimethylnitramine  decomposition  at  pressures  in  the  millitorr  range. 
More  recently,  Lin  and  co-workers  have  studied  the  decomposition  of  dimethyl- 
nitramine  in  both  low  temperature  (466  to  524  K)  static  bulb  experiments  and 
in  single  pulse  shock  tube  work  at  much  higher  temperatures  (800  to  1200  K).^ 

C 

Lee  and  coworkers  have  also  recently  studied  the  infrared  multiphoton  decompo¬ 
sition  of  dimethylnitramine  under  collisionless  conditions.  The  results  of 
these  and  other  workers  are  summarized  in  Table  1.  In  all  but  the  two  cases 
where  the  molecules  are  the  most  isolated  (the  multiphoton  and  the  VLPP 
studies),  the  major  product  is  dimethylnitrosamine .  Thus,  reaction  leads  to  a 
product  having  the  N-N  bond  intact,  in  spite  of  the  fact  that  the  initial  reac¬ 
tions  are  presumed  to  have  involved  N-N  bond  scission  of  one  sort  or  another. 
Reforming  the  N-N  bond  (e.g.,  through  radical  recombination  with  NO)  is  not  in 
itself  difficult,  but  for  it  and  other  secondary  reactions  to  take  place  at 


rates  such  that  dimethylnitrosamine  is  produced  in  almost  quantitative  yield  is 
a  result  that  one  would  not  have  predicted. 
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The  general  goal  of  this  work  is  to  determine  the  products  of  dimethylni- 


tramine  decomposition  and  the  temperature  dependence  of  the  various  rate- 
controlling  processes,  in  order  to  help  achieve  the  overall  goal  of  understanding 
what  controls  the  decomposition  of  cyclic  nitramines  and  other  energetic 
materials.  Under  almost  all  circumstances  (all  but  low  concentrations  and/or 
very  short  time  scales) ,  the  nature  of  these  materials  dictates  that  this  means 
considering  not  only  the  initial  steps  but  also  a  varying  range  of  secondary 
reactions.  This  paper  presents  our  results  for  the  thermal  decomposition  of 
dimethylnitramine  and  dimethynitrosamine  using  GC-MS  monitored  laser-powered 
homogeneous  pyrolysis  (LPHP) .  In  this  technique,  the  net  decomposition  in  a 
flow  system  averaged  over  many  laser  pulses  and  over  the  entire  volume  of  the 
cell  is  measured.  A  mechanistic  numerical  model  accounting  for  the  observed 
Arrhenius  parameters,  the  products,  and  their  relative  yields  is  also 
presented.  Complementary  results  obtained  with  a  real-time  molecular  beam  mass 
spectrometrically  monitored  version  of  LPHP  are  the  subject  of  a  paper  to  be 
submitted. 


EXPERIMENTAL  PROCEDURE 


The  pulsed  laser  pyrolysis  technique  used  here  (GC/MS-monitored)  has  been 

described  in  detail.  Briefly,  a  CO2  laser  is  used  to  heat  the  substrate  (at 

0.1  torr)  via  an  absorbing  but  unreactive  gas  (e.g.,  SFg  -10  torr) ,  which  then 

transfers  its  kinetic  energy  by  collision  (total  pressure  100  torr,  collision 

q 

frequency  -  10  /s) ,  bringing  the  substrate  to  high  temperature  in  1-3  ns.  For 
radical  forming  processes,  scavengers  and/or  traps  are  necessary  (e.g.  toluene, 
cyclopentane,  NO  -5-20  torr)  Reaction  then  occurs  in  the  roughly  10  ns 
available  before  cooling  takes  place  as  the  acoustic  expansion  wave  moves 
radially  inward  toward  the  center  of  the  cell.  Figure  1  shows  a  schematic  of 
the  LPHP  apparatus . 

With  the  LPHP  technique,  surface  catalyzed  reactions  at  high  temperatures 

are  essentially  precluded,  since  the  time  scale  for  heating,  reaction,  and 

cooling  is  far  shorter  than  the  millisecond  time  scale  of  diffusion  to  the 

walls.  (Adsorption  of  IR-absorbing  reactants  or  products  on  the  KC1  windows 

could  conceivably  result  in  decomposition  of  an  adsorbed  phase.  This  is 

generally  unlikely,  since  only  extremely  strongly  absorbing  materials  would  be 

heated  to  a  sufficiently  high  temperature  when  present  at  only  a  few  monolayers 

thickness  on  the  transparent,  cold  windows.) 

The  need  to  explicitly  measure  the  temperature  corresponding  to  any 

particular  measurement  of  k  is  eliminated  by  use  of  the  "comparative  rate" 

8  9 

technique.  ’  This  involves  concurrent  determination  of  the  fractional 
decomposition  of  a  temperature  standard,  a  second  "substrate"  (-0.2  torr)  whose 


decomposit ‘ ->n  rate  temperature  dependence  is  already  well  known,  and  therefore 
whose  fractional  decomposition  defines  an  "effective"  temperature  that  is  a 
convolution  of  temperature  over  time  and  space  within  the  heated  portion  of  the 
cell.  Given  that  the  temperature  dependence  of  the  substrate  and  the  tempera¬ 
ture  standard  are  not  greatly  different,  and  that  the  fractional  decomposition 
of  both  materials  is  less  than  about  25%,  the  standard  and  the  substrate 
"track"  the  temperature  variations  in  a  similar  manner.  Data  reduction  simply 
involves  determining  the  slope  and  intercept  of  a  plot  of  log  kunknown  vs  l°g 
^standard’  correspond  to  the  ratio  of  activation  energies  and  the  differ¬ 

ence  in  the  log  A  values,  respectively.  As  discussed  in  references  8  and  9, 
uncertainties  in  reaction  time  and  in  cell  volumes  cause  no  error  in  the 
derived  activation  energies  (slope  of  the  comparative  rate  plot)  and  only  small 
errors  in  derived  A- factors. 

The  extent  to  which  the  comparative  rate  method  accommodates  the  tempera¬ 
ture  and  reaction  time  variations  inherent  in  this  laser  pyrolysis  technique 
has  been  amply  demonstrated  by  estimates  of  the  errors  involved  and  by  testing 
the  technique  using  a  series  of  "unknowns"  whose  Arrhenius  parameters  have 
already  been  reported  in  the  literature.  ’  This  has  provided  values  in  good 
agreement  with  the  literature  parameters  for  a  series  of  molecular  elimination 
reactions.  The  precision  achievable  indicates  that  the  pulse-to-pulse  and 
minute -to -minute  laser  variations  that  existed  under  the  conditions  of  this 
experiment  did  not  give  rise  to  any  distortion  of  the  derived  temperature 
dependence.  Tests  of  this  type  make  it  clear  that  the  comparative  rate  method, 
applied  to  pulsed  laser  pyrolysis,  is  very  successful  in  accommodating  *"he 
temperature  variations  inherent  in  the  method  and  thereby  providing  an  accurate 
measure  of  the  temperature  dependence  of  the  substrate. 


In  this  particular  study,  initial  experimental  difficulties  included 
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adsorption  of  substrates  on  the  cell  walls,  premature  decomposition  of 


substrates  and  interference  by  chemical  impurities.  These  difficulties  were 


minimized  and  good  mass  balances  obtained  by  silanizing  the  entire  flow  system 


including  the  gas  reservoirs,  carefully  maintaining  the  lines  and  cell  tempera¬ 


tures  at  95°C  (at  lower  tempertaures ,  adsorption  occurred;  at  higher  tempera¬ 


tures,  preliminary  decomposition  occurred),  and  by  scrupulous  purification  of 


all  compounds  used.  In  particular,  traces  of  NO2  were  removed  from  NO  by 


passing  it  through  3  scrubbers  containing  10  M  NaOH,  NaOH  pellets,  and 


molecular  sieves,  respectively.  Special  care  was  given  to  the  elimination  of 


NO2  from  NO,  since  the  former  readily  gave  rise  to  the  nitrates  and  nitrites  of 


the  temperature  standards  used,  thus  distorting  their  decomposition  temperature 


dependence , 
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RESULTS 


A.  Experiments 


The  thermal  decomposition  of  dimethylnitramine  was  studied  over  the  tem¬ 


perature  range  780  K  to  880  K.  Figure  2  shows  the  decomposition  rate  of 


dimethylnitramine  versus  that  of  4-hydroxy-4-methyl-2-pentanone  (HMP)  for  two 


different  experimental  conditions.  The  rate  expression10  for  HMP  is 


log  k(HMP  -»  2  acetone)  -  11.63  -32.3/2.3RT.  In  most  cases  an  approximately 


50-  to  100-fold  excess  of  radical  scavenger,  i.e.  either  toluene  or  cyclopen¬ 


tane,  was  used  and  the  mass  balance  (decomposed  material  accounted  for  as 


product)  for  the  temperature  standard  was  99±3%. 


Excess  nitric  oxide  was  used  in  two  of  the  runs  shown  in  Figure  2  in  order 


to  trap  any  dimethylaminyl  radicals  formed  by  the  homolysis  of  the  N-N  bond  in 


dimethylnitramine.  Secondary  reactions  and  products  were  thus  minimized  and 


dimethylnitrosamine  accounted  for  100±16%  of  the  decomposed  dimethylnitramine 


(for  decompositions  that  were  typically  less  than  10%).  The  least  squares  line 


drawn  through  the  filled  and  un- filled  squares  is  a  composite  of  two  runs 


performed  under  similar  conditions.  Toluene  was  used  as  the  radical  scavenger 


with  a  50- fold  excess  of  nitric  oxide  in  one  run  (filled  squares)  and  a 


100-fold  excess  in  the  other  (unfilled  squares).  As  shown  in  Table  2,  the 


Arrhenius  parameters  for  the  dimethylnitramine  decomposition  in  the  presence  of 


excess  NO,  using  either  toluene  or  cyclopentane  as  the  radical  scavenger,  are 


essentially  the  same,  i.e.,  log  k(s*i)  -  13.5  ±  0.6-(37.4  ±  2.5)/2.3RT. 
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Data  obtained  with  toluene  as  the  scavenger  and  in  the  absence  of  nitric 
oxide  are  also  plotted  in  Figure  2.  In  this  case  no  more  than  60%  of  the 
decomposed  dimethylnitramine  could  be  accounted  for,  appearing  as  diraethylni- 
trosamine  (25%),  methylmethyleneimine  (-20%),  nitrous  acid  (-12%),  dimethyl- 
amine  (-8%),  and  N, N- dime thy lbenzylamine  (-3).  These  yields  varied  with 
temperature,  with  the  dimethylnitrosamine  yield  being  smaller  at  the  lower  end 
of  the  temperature  range.  Trace  amounts  of  N,N-dimethylacetamide  were  tenta¬ 
tively  identified.  Unlike  runs  with  excess  NO,  in  which  virtually  identical 
comparative  rate  plots  are  obtained  either  from  reactant  disappearance  or 
product  appearance,  the  comparative  rate  plot  for  data  obtained  in  the  absence 
of  excess  NO  could  only  be  based  on  the  amount  of  dimethylnitramine  lost. 
Nevetheless,  as  shown  in  the  figure  and  in  Table  2,  the  Arrhenius  parameters 
obtained  for  dimethylnitramine  decomposition  in  the  absence  of  excess  nitric 
oxide  are  still  similar  to  those  obtained  in  its  presence. 

Prior  to  using  HMP  as  the  temperature  standard  in  the  thermal  decomposi¬ 
tion  of  dimethylnitramine,  t-butyl  bromide  was  used  as  the  temperature 
standard^  in  several  runs  In  which  the  nature  of  the  radical  scavenger  or  of 
the  nitrogen  oxide  present  in  the  gas  mixture  was  varied.  The  Arrhenius  param¬ 
eters  obtained  from  the  resulting  comparative  rate  plots  were  similar  to  those 
obtained  with  HMP.  However,  due  to  experimental  difficulties  (see  experimental 
section) ,  attempted  mass  balances  for  both  the  temperature  standard  and  the 
dimethylnitramine  ranged  widely.  The  quantitative  results  were  therefore 
considered  unreliable  and  are  not  reported  here  (even  though  the  Arrhenius 
parameters  generally  fall  within  the  range  shown  in  Table  2).  The  following 
qualitative  observations  from  these  runs  are  nonetheless  considered  pertinent 


in  understanding  this  system's  complexities  and  we  attempted  to  incorporate 
them  in  the  mechanistic  modeling  of  the  DMNA  decomposition. 

In  the  absence  of  excess  nitrogen  oxides,  use  of  toluene  as  a  radical 
scavenger  resulted  in  the  formation  of  N.N-dimethylbenzylamine ,  benzene  and 
bibenzyl  in  addition  to  the  formation  of  N,N-dimethylnitrosamine  as  the  major 
product  and  formaldehyde  in  trace  amounts.  Use  of  cyclopentane  as  the  radical 
scavenger  gave,  in  addition  to  nitrosamine  and  formaldehyde,  nitrocyclopentane 
and  cyclopentene . 

Since  the  reduction  of  NO2  in  secondary  reactions  is  considered  to  be  a 
critical  step  in  dimethylnitrosamine  production  via  recombination  of  the 
dime  thy lamino  radical  with  NO,  dimethylnitramine  decomposition  was  also  studied 
in  the  presence  of  toluene,  t-butylbromide  and  roughly  a  100- fold  excess  of 
nitrogen  dioxide.  Although  a  temperature  standard  was  not  included  in  this 
run,  excess  NO2  appeared  to  have  no  substantial  effect  on  the  dimethylnitramine 
decomposition  rate. 

When  the  dimethylnitramine  decomposition  was  carried  out  in  the  absence  of 
radical  scavenger  or  of  nitrogen  oxides,  dimethylnitrosamine  continued  to  be 
the  major  product,  but  formaldehyde  was  produced  in  more  than  50%  yield 
(relative  to  nitrosamine  produced),  and  the  rate  of  nitramine  loss  apparently 
increased  by  a  factor  of  from  2  to  5.  (The  exact  extent  of  acceleration  is 
uncertain  because  the  temperature  standard  was  also  absent  in  these  runs.) 

In  order  to  check  the  behavior  of  our  temperature  standard,  HMP,  its 
Arrhenius  parameters  were  independently  determined  using  the  LPHP  comparative 
rate  technique  with  t-butylbromide  as  the  temperature  standard  (literature 
parameters11  log  k/s'1  -  13.8  -  41.8/2.3RT).  The  A-factor  and  activation 
energy  for  HMP  thus  determined  were  in  excellent  agreement  with  those  reported 


in  the  literature^-®  (log  A  -  11.9  vs  11.6,  and  Ea  -  32.9  vs  32.8  kcal/mol, 
respectively).  This  result,  along  with  the  good  mass  balances  obtained  for 
diraethylnitramine  decomposition  in  the  presence  of  excess  nitric  oxide  and  the 
similar  Arrhenius  parameters  found  for  dimethylnitramine  using  t-butyl  bromide 
as  the  temperature  standard,  all  suggest  that  the  Arrhenius  parameters  we 
report  for  dimethylnitramine  are  not  distorted  by  systematic  errors  resulting 
from  different  temperature  tracking  behaviors. 

Because  of  its  importance  as  a  product,  the  thermal  decomposition  of 
dimethylnitrosamine  was  also  studied  over  the  temperature  range  870  to  990  K, 
using  isopropyl  bromide  as  the  temperature  standard.  The  comparative  rate  plot 
is  shown  in  Figure  3.  With  dimethylnitrosamine  the  sequence  of  secondary  gas 
phase  reactions  is  considerably  simpler  than  with  dimethylnitramine  since  the 
formation  of  imine  by  loss  of  a  hydrogen  atom  is,  in  the  presence  of  scavenger, 
essentially  irreversible,  and  the  sequence  of  oxidation  steps  that  take  place 
when  NO2  is  present  cannot  occur  here.  However,  GC  analysis  of  the  imine 
presents  a  problem,  since  the  proclivity  of  imines  for  polymerization  reactions 
means  that  any  surface  contacted  by  the  imine  (e.g.,  in  the  heated  transfer 
lines)  is  a  candidate  for  promoting  polymerization.  The  extent  of  this  poly¬ 
merization  varies  from  day  to  day,  exhibiting  a  sens  ivity  to  factors  not 
easily  defined.  Samples  of  the  gummy  polymer  were  collected  and  subjected  to 
analysis  by  field  ionization  mass  spectrometry  (FIMS).  The  dominant  peak  in 
the  pyrolysis/FIMS  spectra  of  this  material  was  m/z  86,  corresponding  to  an 
imine  dimer  produced  from  thermal  decomposition  of  the  imine  polymer. 

The  net  result  of  this  behavior  is  that,  although  the  gas  phase  secondary 
reactions  of  dimethylnitrosamine  are  much  simpler  than  for  dimethylnitramine, 
reliable  quantitative  product  analysis  is  very  difficult.  Hence  the  kinetics 
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were  based  solely  on  DMNO  disappearance,  and  there  is  substantial  scatter  in 

the  comparative  rate  plot.  This  scatter  notwithstanding,  the  parameters 

extracted  from  the  data  in  Figure  3  are  consistent  both  with  the  expectation  of 

simple  N-NO  bond  scission  as  the  rate-limiting  step,  and  with  the  BAC-4  calcu- 

lations  of  Melius  that  suggest  N-NO  bonds  are  not  weaker  than  their 

N-NO2  counterparts,  in  marked  contrast  to  C-NO  bonds  and  their  C-NO2 

counterparts.  Thus,  this  result  serves  to  reiterate  that  the  LPHP  technique 

behaves  as  claimed,  and  impels  us  to  examine  the  chemistry  more  closely  in  the 

case  of  dimethylnitramine.  J 


B.  Numerical  Modeling  I 

In  considering  the  data  presented  above,  the  goal  was  to  choose  a  sequence  ! 

of  reactions  that:  (1)  explained  the  observed  products  and  the  relative  yield  1 

| 

in  which  they  are  formed,  (2)  explained  the  observed  overall  rate  parameters,  \ 

i 

and  (3)  did  so  with  known  or  plausible  parameters  for  each  of  the  elementary  j 

processes.  The  mechanistic  numerical  model  we  constructed  to  help  explain  the 
results  (decomposition  of  dimethylnitramine)  consists  of  75  primary  and  secon¬ 
dary  reactions  that  account  for  all  of  the  observed  products.  Reverse  reac¬ 
tions  were  initially  included  for  all  reactions;  some  were  dropped  from  the 
model  when  thermochemically  limited  parameters  and  numerical  Integration  showed 
them  to  be  of  negligible  importance.  Certain  reactions  postulated  in  the 
literature  as  important  contributors  for  the  dimethylnitramine  decomposition 
were  retained  in  the  model  even  if  thermochemical  limitations  ruled  out  their 
importance . 

J  The  Arrhenius  parameters  for  all  reactions  were  taken  from  the  literature 

1 


where  possible;  otherwise  they  were  derived  by  estimation  techniques  from  anal¬ 
ogous  reactions  in  the  literature.  Where  appropriate,  the  degree  of  falloff 


was  estimated  by  the  RRK  method.  Instead  of  discussing  each  reaction  and  the 
parameters  used  for  it  individually,  the  main  reactions  in  the  model  are  sche¬ 
matically  depicted  in  Figure  4.  In  Table  3  we  list  the  rate  parameters  for  the 
reactions  included  in  Figure  4.  We  have  endeavored  to  assure  maximum  utility 
for  the  modeling  effort  by  remaining  cognizant  of  the  fact  that  literature 
data^  for  all  of  the  reaction  types  allow  only  a  limited  range  within  which 
the  parameters  can  be  varied  without  pressing  the  limits  of  plausibility. 

In  order  to  approximate  the  physical  characteristics  of  the  LPHP  system 
and  to  account  for  the  possibility  of  reaction  during  cooling  (primarily  low 
activation  energy  secondary  reactions  of  radicals)  of  the  laser-heated  portion 
of  the  cell,  the  model  incorporated  a  drop  of  200K  (from  the  initial  900K)  at 
the  end  of  the  first  time  decade  (10  ps) ,  followed  by  100K  drops  at  the  end  of 
each  succeeding  decade .  The  sequence  of  differential  equations  provided  by  the 
reactions  was  then  integrated  over  time  and  temperature  by  a  program  that  uses 
a  numerical  integration  routine  based  on  the  Gear  algorithm.  The  temperature 
profile  and  the  computed  dimethylnitramine  concentration  as  a  function  of  time 
is  shown  in  Figure  5. 

The  experimental  temperature  dependence  and  product  yields  under  different 
conditions  were  successfully  reproduced  only  when  the  model  included  the  nitro- 
nitrite  rearrangement  and  also  reactions  4,  6,  and  11.  These  reactions  are 
reversal  of  the  nitro-ni trite  rearrangement  and  of  the  nitrite  thermolysis,  and 
displacement  of  NO2  from  the  nitrite  by  NO,  respectively.  The  reasonableness 
of  these  reactions  is  considered  in  the  following  section;  taken  at  face  value, 
the  computed  values  suggest  that  the  experimental  results  do  in  fact  reflect 
initial  competition  between  simple  bond  scission  and  nitro-nitrite  rearrange¬ 
ment,  each  with  the  parameters  shown  in  Table  3. 


DISCUSSION 


As  described  in  the  preceding  section,  the  mechanism  whose  important 
elements  are  shown  in  Figure  4  is  able  to  reproduce  the  data,  including  nitra- 
mine  decomposition  rate  parameters  and  dimethylnitrosamine  yield  under  a  range 
of  experimental  conditions,  and  furthermore  does  so  with  rate  parameters  for 
the  various  steps  that  are  either  experimentally  determined  or  plausible  extra¬ 
polations  from  related  reactions.  In  this  section,  we  examine  the  key  steps  in 
the  model  and  discuss  them  relative  to  those  that  have  been  suggested  in  prev¬ 
ious  studies  to  account  for  observed  rates  and  products. 

The  basic  experimental  results  that  must  be  duplicated  if  a  dimethylnitra- 
mine  decomposition  mechanism  is  to  be  considered  potentially  valid  are  the 
decomposition  rate  and  dimethylnitrosamine  yield  as  a  function  of  temperature 
and  radical  scavenger/trap  concentrations.  The  difficulty  that  this  seemingly 
simple  requirement  presents  is  that  the  experimental  parameters  are  too  low  for 
simple  bond  scission  to  be  the  dominant  rate -determining  step,  but  any  molecu¬ 
lar  elimination  or  rearrangement  pathway  that  one  writes  to  allow  duplication 
of  the  Arrhenius  parameters  does  not  lead,  by  a  simple  sequence  of  known  reac¬ 
tions,  to  the  observed  dimethylnitrosamine  yields.  Specifically,  such  sequen¬ 
ces  do  not  provide  a  route  from  the  nitroxyl  radical  (which  is  the  unavoidable 
scission  product  that  follows  nitro-nitrite  rearrangement)  to  the  nitrosamine.* 


One  might  have  imagined  that  HONO  from  molecular  elimination  could  take  part 
in  some  surface-promoted  nitrosation  of  amine  during  transport  of  the  reaction 
products  to  the  GC/MS  system,  but  this  would  have  involved  reaction  with  the 
amine,  rather  than  the  imine.  Furthermore,  the  imine  tends  to  undergo  a 
variable  amount  of  polymerization  on  the  walls  of  the  the  tranfer  lines;  the 
observed  yield  of  dimethylnitrosamine  is  quite  reproducible,  and  therefore  not 
easily  connected  with  reactions  subject  to  variations  in  surface  conditions. 


Given  the  notorious  susceptibility  of  even  precisely  measured  Arrhenius 
parameters  to  distortion  by  systematic  error,  we  were  initially  inclined  to 
assign  the  low  temperature  dependence  to  such  distortions  and  conclude  that  the 
100%  yield  of  nitrosamine  produced  in  the  presence  of  excess  NO  simply  resulted 
from  trapping  a  quantitative  yield  of  dimethylamino  radicals  produced  in  an 
initial  N-NO2  bond  scission.  However,  rate  measurements  made  under  a  wide 
range  of  experimental  conditions  that  were  employed  in  an  effort  to  minimize 
distortion  of  the  temperature  dependence  by  spurious  secondary  reactions  have 
repeatedly  given  log  k(s*^)  -  13.5  ±  0.5  -  (37.4  ±  1.5)/2.3RT,  as  described  in 
the  preceding  section.  Furthermore,  we  conclude  that  there  is  no  significant 
distortion  of  the  temperature  dependence  by  the  highly  reactive  imine,  since 
pyrolysis  of  dimethylnitrosamine ,  which  yields  the  imine  as  the  dominant 
product,  produces  a  temperature  dependence  completely  consistent  with  an  N-NO 
bond  scission.  [Finally,  preliminary  results  of  dime thy Ini tramine  pyrolysis  in 
a  molecular -beam  mass  spectrometrically  sampled  LPHP  cell  reveal  that  mass  60 
and  mass  30  (presumably  dime thy Ini troxyl  radical  and  NO)  are  products  of  an 
initial  unimolecular  decomposition  step.]  Thus,  the  results  indicating  that  a 
low  A- factor  initial  step  contributes  to  dimethylnitramine  decompostion  became 
too  consistent  to  simply  reject  out  of  hand,  and  demanded  a  plausible  mechanism 
that  could  account  for  them. 

The  critical  reactions  that  make  it  possible  for  the  numerical  model  to 
duplicate  the  low  Arrhenius  parameters  and  high  dimethylnitrosamine  yields  in 
the  presence  of  excess  NO  are  shown  in  Figure  4  and  Table  3,  which  have  been 
abridged  for  purposes  of  discussion.  The  critical  elements  are: 

1.  N-NO9  bond  scission  and  NO2-ONO  rearrangement  that  are 
roughly  competitive  at  900  K  (Reactions  1  and  3). 


2.  Rapid,  but  reversible  thermolysis  of  the  very  weak  0-N0 
bond  in  the  N-nitrite  (Reactions  5  and  6) . 

3.  Reverse  of  the  nitro-nitrite  rearrangement  (Reaction  4). 

4.  Attack  of  NO  on  the  N-nitrite  to  displace  NO2  and  yield 
directly  the  nitrosamine  (Reaction  11) . 

The  reverse  of  the  nitro-nitrite  rearrangement  and  of  the  nitrite  thermolysis 
are  required  by  microscopic  reversibility  (though  we  did  not  initially  antici¬ 
pate  their  impact).  Reversal  of  the  0-N0  bond  thermolysis  is  important  only 
in  the  presence  of  excess  NO,  and  it  is  only  this  reversal  that  can  maintain  a 
high  enough  steady  state  concentration  of  the  weakly  bonded  nitrite  such  that 
reversal  of  the  nitro  nitrite  rearrangement  (Reaction  4)  becomes  significant. 

It  is  this  latter  reaction  that  allows  some  return  to  the  nitramine  and  helps 
to  limit  the  formation,  under  these  conditions,  of  other  products  from  the 
nitroxyl  radical. 

The  only  reaction  in  the  sequence  without  direct  precedent  is  reaction  (11), 
the  displacement  of  NO2  from  the  nitrite  by  NO.  This  reaction  is  most  impor¬ 
tant  at  lower  temperatures.  Because  it  is  a  low  activation  energy  process  and 
it  limits  return  from  the  nitrite,  it  serves  to  substantially  increase  the  net 
decomposition  at  low  temperatures,  while  providing  an  additional  route  to 
nitrosamine.  As  an  addition-elimination  reaction  between  a  free  radical  and  a 
substrate  having  partial  double  bond  character,  it  is  the  type  of  reaction  that 
could  well  proceed  with  a  low  activation  energy.  As  currently  incorporated  in 
the  model,  it  has  a  5  kcal/mol  activation  energy.  If  this  is  raised  to  10 
kcal/mol,  the  computed  overall  activation  energy  for  DMNA  loss  is  raised  from 
39.5  to  44.5  kcal/mol  (i.e.,  75%  of  the  way  to  the  46.1  kcal/mol  that  results 
from  complete  elimination  of  reaction  11).  There  may  conceivably  be  alterna¬ 
tives  to  reaction  (11);  as  the  reaction  sequence  is  currently  written,  this 
reaction  constitutes  a  necessary  and  reasonable  route  to  nitrosamine  from  the 


easily  formed  nitroxyl  radical.  That  is,  if  nitrosamine  is  to  be  formed  from 


the  nitroxyl  radical,  the  oxygen  on  the  nitrogen  clearly  must  depart,  and 


departing  as  one  of  the  oxygens  on  an  NO2  seems  an  emminently  reasonable  way  to 


accomplish  this.  It  is  interesting  to  note  that  the  stability  of  the  nitroxyl 


radical  at  the  same  time  makes  it  the  inevitable  product  of  the  nitro-nitrite 


rearrangement  and  a  bottleneck  in  the  further  decomposition  of  the  amine  frame¬ 


work.  This  bottleneck,  in  turn,  is  what  evidently  makes  it  possible  for  a 


major  ultimate  product  of  the  nitro-nitrite  rearrangement  to  be  the  nitrosamine. 


In  this  and  the  following  paragraphs,  we  consider  some  of  the  other  reac¬ 


tions  that  have  been  previously  suggested  for  nitrosamine  formation.  Beginning 


with  the  results  of  Fluornoy^  some  years  ago,  a  number  of  suggestions  have  been 


made  for  the  formation  of  dimethylnitrosamine  as  the  main  product  from  the 


thermal  decomposition  of  dimethylnitramine .  Two  of  these,  namely  unimolecular 


oxygen  atom  loss  (15),  and  NO  addition  to  the  imine  followed  by  scavenging  of 


the  resulting  carbon-centered  radical  (16),  can  be  ruled  out  on  thermochemical 


grounds . 
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Reaction  (15)  is  not  a  possible  explanation  because  it  has  an  endothermicity 


substantially  exceeding  the  observed  activation  energy  for  nitrosamine 


production. 
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Reaction  (16/16a)  is  unfavorable  because  the  exothermic  decomposition  of  the 


intermediate  Reaction  (-16)  is  much  faster  than  its  stabilization  by  any 


1 

conceivableiJ  scavenger. 


The  remaining  candidates  that  need  to  be  considered  include  recombination 


of  NO  with  the  dimethylamino  radical  (9)  as  invoked  by  Lin,  radical  addition 


to  one  of  the  nitro-group  oxygens  as  discussed  by  Melius1^  (12,13)  (i.e., 


direct  biotnolecular  reduction  of  the  intact  nitramine),  and  bimolecular  or 


unimolecular  formation  of  the  unstable  O-nitroso  compound  (6,3)  as  an  effective 


route  to  high  yields  of  NO  to  serve  as  an  intermediate  for  reaction  (9). 
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The  recombination  of  NO  with  the  dimethylamino  radical  is  the  logical 
first  candidate  and  has  been  invoked  by  all  previous  workers  as  a  major  source 
of  nitrosamine;  this  step  can  be  preceded  by  NO  production  routes  consisting  of 
a  sequence  of  "known"  unimolecular  and  bimolecular  reactions  (1  and  7)  and/or 
by  the  nitro-nitrite  rearrangement  (3)  which  has  precendent  as  a  unimolecular 
process  in  ni trome thane ^  and  nitrobenzene^ ’ ^  decomposition. 

The  four-center  "recombination"  leading  directly  to  two  nitrosomethane 
molecules  suggested  by  Lin  would  appear  to  have  steric  requirements  that  are 
too  stringent  to  be  the  high  A- factor  process  that  is  evidently  required  to 
explain  his  observations . 


On  the  other  hand,  Reaction  (17)  is  not  an  unreasonable  suggestion  on  enthalpic 
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grounds,  given  the  known  tendency  of  C-nitroso  compounds  to  dimerize. 

The  model  used  by  Lin  successfully  reproduced  his  observed  DMNA  decomposi¬ 
tion  and  DMNO  formation  rates  through  a  route  involving  almost  complete  oxida¬ 
tive  fragmentation  of  the  amine  structure  and  formation  of  NO  as  an  intermed¬ 
iate  in  rather  large  yields  in  the  process.  However,  there  are  several  key 
steps  whose  parameters  need  to  be  carefully  assessed  in  order  to  judge  the 
adequacy  of  the  reaction  scheme  as  assembled  by  Lin.  These  are  the  four-center 
recombination  (17)  considered  above,  and  the  oxidation  of  the  dimethylamino 
radical  by  NO2  (7),  and  the  thermolysis  of  nitrosomethane. 


Reaction  (17)  can  be  viewed  as  a  disproportionation  with  rather  severe 
steric  requirements.  While  the  the  impact  of  these  requirements  on  the  rate 


1  R 

may  be  tempered  by  the  fact  that  the  nitro some thane  dimer  is  bound  by  about 
17  kcal/mol,  an  A- factor  of  3  X  10^  1  m’*"  s''*'  is  still  surprisingly  high. 
Furthermore,  for  this  reaction  to  take  place  to  the  exclusion  of  the  recombina¬ 
tion  at  an  oxygen  atom  site  (i.e.,  reaction  7),  a  known  reaction  that  has  no 
unusual  steric  requirements  and  leads  to  oxidation  of  the  radical,  is  even  more 
surprising.  In  principle,  this  "known"  route  can  lead  to  the  same  products  as 
the  four-center  recombination  (17),  since  oxidation  of  the  dimethyl  amino 
radical  gives  a  nitroxyl  radical  (7  ),  which,  upon  ^-scission,  gives  methyl 
radical  and  a  molecule  of  nitrosomethane . 

(CH3)2N-0-  - >CH3N0  +  CH3‘  (18) 

Thermolysis  of  nitrosomethane  then  gives  NO  and  another  methyl  radical,  which 
can  go  on  to  reduce  more  N02  to  NO. 

CH3-NO  - ►  CH3’  +  NO  (19) 

Whether  the  sequence  of  reactions  (7),  (18),  and  (19)  will  in  practice  give  the 

same  results  as  (17)  and  (19)  is  another  question,  since  nitroxyl  radicals  are 

12 

quite  stable  in  general,  and  this  particular  one  is  estimated  to  have  a 

methyl-N  bond  strength  of  44  kcal/mol. 

Finally,  the  unimolecular  dissociation  (19)  of  the  nitrosomethane  itself 

is  another  step  at  which  a  bottleneck  must  be  avoided,  if  the  recombination 

route  to  nitrosomethane  is  to  satisfactorily  explain  the  various  experimental 

2 

results.  The  parameters  that  Lin  has  evidently  found  necessary  for  this  step 
to  be  fast  enough  include  an  A- factor  of  5  X  10*"^  1  m"*-  s'*-.  Values  in  the 
literature*-^  for  loss  of  diatomic  fragments,  as  well  as  consideration  of  the 


overall  entropy  change  in  this  particular  reaction,  indicate  that  5  X  10^  is 
about  an  order  of  magnitude  too  large  for  nitrosomethanc  thermolysis,  without 
consideration  of  the  effects  of  falloff.  In  fact,  use  of  all  the  other 
parameters  chosen  by  Lin,  but  with  the  A- factor  for  this  reaction  lowered  to 
5  X  1015,  lowered  the  computed  dimethylnitrosamine  yield  at  500  K  from  the 
observed  60%  to  40%.  These  difficulties  do  not  absolutely  rule  out  NO 
recombination  as  the  sole  route  to  the  nitrosamine,  but  serve  at  least  to 
(a)  quantify  the  intuition  that  a  delicate  balance  among  secondary  reactions  is 
required  in  order  to  produce  enough  NO  for  nitrosamine  production  through 
recombination,  while  still  leaving  enough  intact  dimethylamino  radicals  to 
recombine  with  the  NO,  and  (b)  amplify  Lin's  original  finding  that  a  fully 
satisfactory  explanation  does  not  easily  follow  from  previously  documented 
reactions . 

A  nitrosamine  formation  route  that  would  not  require  a  critical  balance 

between  amino  radical  survival  and  oxidative  fragmentation  is  the  one  first 
19  12 

made  by  Schroeder  and  recently  subjected  by  Melius  to  calculation  by  the 
BAC-MP4  method.  If  Reaction  (12)  is  facile  enough,  it  will  result  in  rapid 
formation  of  the  nitrosamine  from  virtually  any  radical  produced  in  the  system. 

Reactions  12  (and  13)  are  quite  reasonable:  they  are  molecular  analogs 
to  the  oxidation  of  radicals  by  NO2.  Since  they  are  not  radical -radical 
reactions,  but  radical  additions  to  a  double  bond,  one  might  expect  them 
to  exhibit  activation  energies  typical  of  those  classes  of  reactions 
(1-2  keal/mol  for  H-atoms,  and  5-7  for  carbon-centered  radicals).  However,  the 
calculations  of  Melius  evidently  suggest  little  or  no  activation  energy  for 
methyl  radicals  as  well  as  for  hydrogen  atoms.  Thus,  in  order  to  derive  an 
upper  limit  to  the  possible  importance  of  these  reactions,  we  have  incorporated 
them  (for  R  -  H‘  and  CHj')  into  the  numerical  model  described  above  using 


appropriate  A- factors  and  zero  activation  energies  for  methyl  radicals  as  well 
as  hydrogen  atoms.  Under  these  conditions,  the  computation  showed  that 
reaction  under  the  modeled  temperature -time  history  illustrated  in  Figure  5  and 
in  the  presence  of  a  50-fold  excess  of  toluene  as  radical  scavenger  resulted  in 
a  contribution  to  diraethylnitrosamine  production  from  reactions  12  and  13  that 
was  no  more  than  20%  of  the  total  yield  of  that  product.  Furthermore,  the 
inclusion  of  these  reactions  did  not  lead  to  a  tendency  towards  lower  values. 
Thus,  while  they  may  be  important  during  actual  nitramine  combustion,  the 
nitrosamine  formation  reactions  considered  by  Melius  would  seem  not  to  be 
responsible  for  the  low  A-factors  and  low  activation  energies  observed  under 
our  scavenged  laser  pyrolysis  conditions.  In  other  words,  we  cannot,  by 
inclusion  of  these  reactions  in  the  model,  duplicate  the  observed  temperature 
dependence  for  DNMA  loss,  namely  an  A- factor  and  activation  energy  two  powers 
of  ten  and  10  kcal/mol ,  respectively,  below  those  indicative  of  simple  N-N02 
bond  scission. 


SUMMARY 


The  decomposition  temperature  dependence  and  the  product  yields  in  the 
laser  pyrolysis  of  dimethylnitramine  can  be  accounted  for  by  a  combination  of 
(1)  competing  unimolecular  bond  scission  and  nitro-nitrite  rearrangement  as  the 
initial  loss  reactions,  and  (2)  by  a  low  temperature  pathway  to  nitrosamine 
that  involves  reversals  of  the  nitro-nitrite  rearrangement  and  nitrite  homoly¬ 
sis  steps,  coupled  with  a  direct  route  from  the  unstable  nitrite  (or  its  oxida¬ 
tion  state  equivalent,  the  nitroxyl  radical)  to  dimethylnitrosamine .  Among 
these,  the  latter  is  the  only  reaction  strictly  without  precedent;  even  it,  as 
a  displacement  of  NO2  from  the  nitrite  by  NO  is  somewhat  analogous  to  radical 
addition  to  a  double  bond.  Parameters  assigned  to  it  on  this  basis  allow  the 
experimental  observations  to  be  reproduced  without  abandoning  plausibility. 

These  results  reveal  competing  initial  reactions  and  facile  secondary 
reactions  in  the  decomposition  of  a  simple  nitramine.  Understanding  the  range 
of  possible  reactions  is  clearly  a  necessary  part  of  understanding  the  condi¬ 
tions  required  for  the  decomposition  of  this  and  more  complex  nitramines  to 
become  thermally  self-sustaining.  It  is  entirely  fitting  that  our  attempt  to 

further  this  understanding  is  presented  as  part  of  a  tribute  to  Fred  Kaufman, 
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who  so  aptly  described  endeavors  such  as  this  in  his  seminal  plenary  lecture 
to  the  19th  Symposium  (International)  on  Combustion.  Experiments  now  underway 
with  a  different  LFHP  apparatus,  in  which  the  initial  products  of  laser 
decomposition  are  detected  directly  via  molecular  beam  mass  spectrometry,  will 
provide  additional  information  about  this  competition. 
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FIGURE  CAPTIONS 


Figure  1  Schematic  of  GCMS  monitored  LPHP  Apparatus. 

Figure  2  Comparative  rate  plot  of  dimethylnitramine  (DMNA)  vs  4-hydroxy-4- 

methyl-2-pentanone  (HMP) .  B*50-fold  excess  of  NO;  Q-100-fold  excess 
of  NO;  x-no  additional  NO;  toluene  as  scavenger. 

Figure  3  Comparative  rate  plot  of  dime thylnitr as amine  (DMNO)  vs  isopropyl 
bromide  (iPrBr).  (The  different  symbols  represent  experiments  on 
different  days . ) 

Figure  4  Schematic  representation  of  the  mechanism  at  DMNA  decomposition  in 
the  presence  of  excess  NO. 

Figure  5  Schematic  representation  of  the  temperature  changes  with  time  imposed 
during  numerical  modeling  to  represent  LPHP  conditions.  Also  shown 
is  the  computed  variation  of  DMNA  with  time.  Initial  conditions  for 
model:  cell  temperature  367° C,  PQ(DMNA)  -  0.11  Torr,  PQ(o- 
fluorotoluene) ,  PQ(N0)  -  5.7  Torr,  PQ(hydroxy-methylpentanone)  - 


DMNA  DECOMPOSITION  MECHANISM 


MODEL  RESULTS 


LOG  k  =  13.7 -39.5/ 


Table  1 

Literature  Parameters  for  Dimetnynitramine  Decomposition 


ACs"1) 

E(kcal/mol) 

Technique 

(T/K)  Range 

P/Torr 

Reference/ 

Year 

a 

Molecular  Beam 
MPD 

Lee/19865 

15.9 

44.1 

Bulb/SPST 

466-1000 

474,3-6 

Lin/19854 

16.510.5 

45.4 

Bulb 

466-524 

475 

Lin/19852 

12.4 

37.0 

VLPP 

560-850K 

-10* 3 

McMillen/19793 

13.7 

38.9 

Bulb 

450-510 

-200 

Korsunski/1987^ 

14.1 

40.8 

Bulb 

450-530 

64-400 

Korsunski/19642 

20.1 

53.0 

Bulb 

440-470 

200-750 

Flournoy/19 6 l^ 

aProduct  mass  spectra  and  flight  times  consistent  with  N-NO2  band  scission. 
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Table  2 

Laser  Pyrolysis  Comparative  Rate  Results  for  Dimethylnitramine 


Temp.  Std. 

Scavenger 

NO  x  Fold 
Excess 

logtA/s'1] 

Ea/kcal/mol 

Nitrosamine 

Yield3 

HMP 

cyclopentane 

100 

13.5 

37.4 

113 

HMP 

toluene 

100,  50 

13.5 

37.4 

99 

HMP 

toluene 

0 

12.8 

35.5 

25 

a(Dime thylnitrosamine)/A(Dimethylnitramine) 


Table  3 

Parameters  Used  in  Numerical  Modeling  of  Dimethylnitramine 
and  Dimethylnitrosamine  Decomposition 


No  ■ _ Reaction _ log  Aa _ E/kcal/mol 


1 

(ch3)2n  -  no2  - 

(CH3)2N  +  no2 

15.3 

46.5 

2 

(ch3)2n  -  no2  - 

(CH3)2N  +  no2 

8.8 

-2.0 

3 

(ch3)2n  -  no2  - 

(CH3)2N-0-N-0 

11.5 

31.0 

4 

(ch3)2n  -  no2  «- 

(CH3)2N-0-N-0 

10.5 

11.1 

5 

(CH3)2N-0-N-0  -► 

(CH3)2N-0>  +  NO 

13.7 

9.0 

6 

(CH3)2N-0-N-0  <- 

(CH3)2N-0-  +  NO 

9.9 

0.0 

7 

(CH3)2N-  +  no2  - 

(CH3)2N-0-  +  NO 

10.0 

0.0 

8 

(CH3)2N-  +  no2  «- 

<CH3)2N-0-  +  NO 

9.9 

13.8 

9 

(CH3)2N-  +  NO  - 

(CH3)2N-NO 

9.8 

-2.0 

10 

(CH3)2N-  +  NO  ♦- 

(CH3)2N-NO 

15.5 

50.3 

11 

(CH3)2N-0-N-0  +  NO  -  (CH3)2N-NO 

+  no2 

9.7 

5.0 

12 

(CH3)2N-N09  +  H- 

-  (CH3)2N-NO  + 

OH- 

10.0 

0.0 

13 

(CH3)2N-N02  +  ch3 

•  -♦  (CH3)2N-NO  + 

ch3o 

9.0 

0.0 

‘Units  of  A  are  and  i(ra-s)'^  for  unimolecular  and  bimolecular  reactions, 
respectively. 
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MOLECULAR- BEAM- SAMPLED  LASER  PYROLYSIS  OF 
DIMETHYLNITRAMINE 
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David  M.  Golden  and  Donald  F.  McMillen 
Department  of  Chemical  Kinetics,  Chemical  Physics  Laboratory 
SRI  International,  Menlo  Park,  CA  94025 

ABSTRACT 

A  pulsed  CO2  laser  was  used  to  heat  a  mixture  of  argon,  SFg ,  and  dimethyl - 
nitramine  (DMNA)  in  a  flow  cell,  such  that  the  periphery  of  the  laser-heated 
region  was  immediately  adjacent  to  a  molecular-beam/quadrupole-mass-spectrometer 
sampling  system.  Detection  of  the  product  mixture  in  the  molecular  beam 
revealed  that  initial  N-NO2  scission  was  accompanied  by  a  comparable  amount  of 
the  nitro-ni trite  rearrangement  (followed  by  rapid  scission  of  the  extremely 
weak  (CH-j^NO-NO  bond)  to  produce  directly  NO  and  the  nitroxyl  radical. 
Production  of  this  NO  by  secondary  bimolecular  reactions  can  be  ruled  out 
because  the  time  available  before  the  product  mixture  was  frozen  by  expansion 
was  too  short,  and  because  the  observed  rise  times  of  the  NO  and  NO2  were 
identical.  Calibration  of  the  detection  system  with  laser-heated  NO  and  NO2  at 
the  reaction  temperature  (900  *  50  °C)  revealed  that  the  rearrangement/bond- 
scission  branching  ratio  is  -  0. 7:1.0.  These  results  reveal  that  previous 
studies  of  dime thy Ini tramine  concluding  that  N-NO2  bond  scission  is  the  sole 
rate-limiting  step  are  in  error,  and  suggest  that  "loose"  nitro-nitrite  rear¬ 
rangements  may  be  a  general  phenomenon  in  N-nitro  compounds  as  they  appear  to 
be  in  C-nitro  compounds. 

£ Postdoctoral  Research  Associate. 

^Present  address:  CIBA-GEIGY,  CH-4333  Munchwillen/AG ,  Schweiz,  Switzerland. 
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INTRODUCTION 


The  initial  decomposition  step(s)  of  nitramines  have  been  of  interest  for 


some  Lime  as  one  of  the  important  keys  to  understanding  what  controls  the  ini¬ 


tiation  and  combustion  of  nitramine  explosives  and  propellants,  respectively . 


For  any  energetic  material  this  information  tends  to  be  difficult  to  obtain, 


since  the  utility  of  such  materials  is  based  on  their  tendency  to  undergo 


extremely  rapid  and  exothermic  decomposition  once  certain  kinetic  bottlenecks 


have  been  passed,  thus  obscuring  the  identity  of  the  bottleneck  reaction  (or 


reactions).  Furthermore,  as  has  been  appreciated  for  some  time,  the  function¬ 


ally  useful  nitramines  have  a  multiplicity  of  candidate  initial  steps  that  can 


be  identified  on  the  basis  of  existing  thermochemical  information;  however, 


this  information  is  not  sufficiently  precise  or  reliable  to  allow  an  a  priori 


distinction  among  the  various  candidates.  Finally,  the  situation  is  further 


aggravated  because  secondary  bimolecular  reactions  of  radicals  with  NO2  are 


known  to  be  extremely  fast  and  can  be  expected  to  complicate  attempts  to  disen¬ 


tangle  the  reaction  sequence. 


For  all  of  the  above  reasons,  dimethylnitramine  (DMNA)  has  been  studied  by 


various  workers  as  a  simple  analog  expected  to  exhibit  some  of  the  reactions 


that  are  important  for  the  cyclic  nitramines.  However,  the  apparent  simplicity 


of  DMNA  notwithstanding,  it  has  been  difficult  to  determine  with  certainty  the 


initial  decomposition  step(s) .  Original  work  by  Fluornoy  and  by  Korsunskii 


and  coworkers^  identified  N-NO2  bond  scission  as  the  initial  step,  although 


this  was  a  matter  of  assignment  rather  than  demonstration,  as  the  temperature 
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dependence  of  the  initial  step  was  obscured  to  an  uncertain  degree  by  secondary 
reactions.  More  recently,  Lin  and  coworkers  reported  low- temperature  static 
bulb  experiments  and  high- temperature  single-pulse  shock  tube  results,  which 
when  coupled  together  give  an  activation  energy,  measured  over  a  wide  tempera¬ 
ture  range  that  is  indicative  of  initial  N-NO2  bond  scission.  Similarly, 

Wodtke  and  Lee^  have  made  preliminary  reports  of  infrared  multiphoton  decompo¬ 
sition  in  a  molecular  beam  also  leading  to  the  conclusion  that  N-NO2  is  the 
exclusive  (i.e.,  >  95%)  initial  reaction.  In  contrast  with  these  results  and 
concurrent  with  the  work  of  Lee,  comparative  rate  GC-analyzed  laser  pyrolysis 
experiments  in  this  laboratory  repeatedly  gave  a  temperature  dependence  for 
DMNA  decomposition  that  is  not  consistent  with  initial  N-NO2  bond  scission.'* 

A  choice  between  those  results  supporting  and  those  contravening  N-NO2 
bond  scission  as  the  sole  rate-limiting  step  cannot  be  made  without  reservations 
The  shock- tube  and  static -bulb  work  were  performed  without  radical  scaveng¬ 
ers,  and  there  is  therefore  no  assurance  that  the  decomposition  rates  are  not 
perturbed  by  secondary  reactions.  In  particular,  for  single-pulse  shock-tube 
kinetic  studies  of  hydrocarbon  decomposition,  there  is  a  general  need  for  scav¬ 
engers  to  inhibit  hydrogen  atom  chain  reactions,  when  those  are  structurally 
possible.  In  the  case  of  the  multiphoton  decomposition  of  DMNA  in  a  molecular 
beam,  there  is,  of  course,  no  doubt  that  bimolecular  secondary  reactions  are 
precluded.  However,  in  the  molecular  beam  of  Wodtke  and  Lee,^  the  sole  diag¬ 
nostic  was  the  kinetic  energy  distribution  of  the  m/z  30  signal,  where  the 
diagnosis  that  none  of  the  mass  30  signal  was  actually  NO  coming  from  a  molecu¬ 
lar  rearrangement  depended  upon  the  inability  to  detect,  the  presence  of  a 
component  with  a  slightly  shorter  flight  time  (than  m/z  46  and  the  m/z  30 
signal  that  arises  from  ionization  and  ion  fragmentation  of  NO2) . 


The  literature  data  outlined  above,  its  limitations  notwithstanding, 
together  with  our  own  previous  experience  with  dimithylnitramine ,  led  us  to 
expect  that  N-NC>2  bond  scission  would  be  the  dominant,  and  perhaps  exclusive, 
reaction  at  900  K.  Thus,  when  we  repeatedly  obtained,  in  GC-monitored  laser 
pyrolysis  experiments,^  a  temperature  dependence  indicative  of  a  contribution 
from  a  low  A-factor  process,  we  responded  by  changing  reaction  conditions  in 
search  of  a  suspected  systematic  error.  When  these  efforts,  which  involved 
varying  the  partial  pressures  and  identities  of  the  radical  scavengers  and  the 
temperature  standard,  and  the  use  of  NO  as  a  radical  trap,  all  still  resulted 
in  a  temperature  dependence  too  low  to  indicate  N-NO2  as  the  sole  rate¬ 
determining  step,  we  concluded  that  direct  observation  of  the  initial  step(s) 
was  indicated.  In  this  paper,  we  describe  the  molecular-beam-sampled,  mass 
spectrometric  detection  of  the  initial  products  from  the  SFg-sensitized  laser 
pyrolysis  of  dime thy Ini tramine . 


EXPERIMENTAL  PROCEDURE 


The  thermal  decomposition  of  dimethylnitramine  was  studied  using  a  laser- 
powered  homogeneous  pyrolysis  system  coupled  with  a  molecular-beam/mass 
spectrometer  detection  system.  A  pulsed  CO2  laser  was  used  to  irradiate  a 
flowing  mixture  of  sulphur  hexafluoride  (SFg)  absorber  molecules,  argon  bath 
gas,  and  DMNA  substrate.  Following  collisional  energy  transfer,  the  heated 
DMNA  thermally  decomposed.  A  molecular  beam  of  the  heated  gas  was  formed  by 
pumping  it  through  a  75-pm  diameter  nozzle.  The  beam  was  then  analysed  using  a 
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quadrupole  mass  spectrometer  (QMS),  which  allowed  the  direct  detection  of  init¬ 
ial  reaction  products. 


Flow  System 

The  reaction  mixture  was  prepared  by  establishing  a  flow  of  a  commercially 
prepared  10%  SFg/argon  mixture  (Matheson)  through  a  glass  tube  containing  the 
DMNA.  The  total  pressure  was  maintained  at  a  constant  level  (30  to  90  Torr) 
using  a  pressure  regulator  and  controlling  the  pump-out  rate  of  the  cell.  The 
"reservoir"  tube  containing  the  DMNA  extended  into  the  heated  flow  cell  through 
an  O-ring-sealed  Cajon  fitting  (Figure  1).  In  this  manner,  controlled  amounts 
of  the  nitramine  vapor  were  introduced  into  the  flow.  The  end  of  the  glass 
reservoir  was  positioned  as  close  to  the  laser-heated  region  as  possible  (to 
minimize  surface  effects) .  The  aluminum  reaction  cell  was  maintained  at  a 
temperature  typically  in  the  range  25-35°C.  The  gas  mixture  was  irradiated 
through  the  KC1  windows  at  2.0  Hz  by  a  pulsed  10. 6 -pm  CO2- laser  beam  of  energy 
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<1.8  J/cm  and  l-/is  duration.  A  mechanical  pump  was  used  to  control  the  gas 
flow  at  a  rate  that  insured  complete  renewal  of  the  reaction  mixture  between 
laser  pulses  (0.5  s).  The  laser-heated  region  was  cylindrical  in  shape,  1.1  cm 
in  diameter  and  2. 5 -cm  long.  The  laser  beam  was  positioned  to  irradiate  the 
reaction  mixture  as  close  as  possible  (<  1mm)  to  the  nozzle,  so  that  the 
expansion  wave  generated  by  the  laser  heating  pushed  the  gas  inside  the  nozzle. 
Thus,  a  few  microseconds  after  the  laser  pulse,  the  mixture  was  frozen  by  the 
adiabatic  expansion,  and  radicals  formed  in  the  early  process  of  decomposition 
were  transported  without  undergoing  further  collision  to  the  QMS  in  less  than  a 
millisecond. 

Beam-Sampling  System 

The  supersonic  divergent  nozzle  consisted  of  a  platinum  disk  (0.80  cm  in 
diameter  and  100  /im  thick)  with  a  75-/io-diameter  hole  at  its  center  (Ernest  F. 
Fullam,  Inc.,  Latham,  NY).  This  disk  was  positioned  at  the  apex  of  a  conical 
support  and  held  firmly  in  place  by  a  screw-cap.  A  skimmer  (electroformed 
nickel,  Beam  Dynamics,  Minneapolis,  MN)  with  a  200-/jm  diameter  at  the  throat 
extracted  the  central  portion  of  the  beam.  The  distance  between  the  nozzle  and 
the  skimmer  entrance  could  be  varied  from  a  few  millimeters  up  to  20  mm  by 
external  adjustment  of  a  micrometer  screw.  A  twin-bladed  knife-edged  chopper 
positioned  at  the  entrance  to  Chamber  2  (containing  the  QMS)  allowed  the  option 
of  chopping  the  beam.  The  chopper  operated  at  «  650  Hz. 


Molecular-Beam  and  Mass  Spectrometer-Detection  System 


diffusion  pumps  with  freon-refrigerated  baffles  (Figure  2).  The  second  chamber 
contained  an  Extranuclear  Laboratories  quadrupole  mass  spectrometer  with  elec¬ 
tron  impact  ionization  and  ion-extraction  perpendicular  to  the  beam.  The  posi¬ 
tion  of  the  QMS  was  adjustable  from  the  outside  by  three  screws  and  a  bellows 
suspension.  The  bulk  of  the  neutral  beam  passed  through  the  differentially 
pumped  ionizer  region  and  impinged  on  a  liquid-nitrogen  cold  trap.  The  cold 
trap  prevented  buildup  of  nitramine  in  the  background  gas  during  operation. 

The  entire  pumping  system  could  be  degassed  by  heating  up  to  200°C.  The  ulti¬ 
mate  pressure  (with  no  beam)  was  on  the  order  of  10* ^  torr  in  each  chamber. 


Data  Aquisition 

The  mass -selected  *ons  were  accelerated  into  an  electron  multiplier.  With 
the  infrared  heating  laser  off,  the  CW  molecular  beam  could  be  chopped,  and  the 
signal  for  a  specific  ion  mass  could  be  phase-sensitively  detected  using  a 
lock-in  amplifier  (ITHACO,  Dynatrac  391).  The  QMS  could  then  be  accurately 
tuned  into  any  particular  mass  by  fine  adjustment  of  the  mass  tuning  to  give 
the  maximum  signal  on  the  lock-in  amplifier.  The  QMS  resolution  was  always 
better  than  0.5  amu  and  was  checked  periodically. 

With  the  QMS  tuned  into  a  specific  mass  and  with  the  chopper  off,  the 
heating  laser  was  pulsed.  The  Infrared  light  passed  through  the  cell,  and 
absorption  by  the  SFg  heated  the  contents  of  the  cell.  Laser  light  was  used  to 
trigger  a  1  MHz  (sampling  rate  1  ps  per  point)  transient  digitizer  controlled 
by  a  DEC  LSI  11/23  microprocessor  to  acquire  the  QMS  signal.  To  improve 
signal - to-noise  ratios,  the  signals  from  2000-3000  laser  shots  were  averaged. 

As  a  check  that  no  product  build-up  was  occurring  during  the  experiment,  the 
laser  was  switched  off  and  the  process  repeated.  In  this  case,  a  pulse 


generator  operating  at  50  Hz  was  used  to  trigger  the  digitizer.  The  entire 


process  was  repeated  for  every  mass  studied.  At  the  end  of  the  experiment,  the 


gain  of  the  ion  detector  was  measured  by  observing  the  variation  of  the  CW 


molecular-beam  signal  with  the  high  voltage.  This  "gain  factor"  was  found  to 


be  mass  independent  and  constant  over  long  periods.  The  nitramine  reservoir 


was  then  removed  and  the  background  levels  in  the  vacuum  chamber  for  all  masses 


of  interest  were  measured. 


Data  Reduction 


Three  operations  were  carried  out  on  the  raw  data  for  each  mass  to  yield 


the  final  results.  These  are  briefly  described  below: 


(1)  The  backgrounds  were  subtracted  from  the  raw  data.  This 
operation  corrected  for  all  signals  not  directly  from  the 
beam,  i.e. ,  background  pressure  in  the  ionization  chamber. 

(2)  The  background- corrected  data  was  then  multiplied  by  an 
appropriate  "gain  factor"  corresponding  to  the  high  voltage 
used  when  collecting  the  signals  from  that  mass.  The  mag¬ 
nitude  of  this  factor  was  determined  by  arbitrarally 
selecting  one -voltage  as  having  a  gain  factor  of  1.0  and 
then  calculating  the  factor  for  all  other  voltages.  This 
procedure  corrected  for  the  signal -enhancement  effect  of 
raising  the  high  voltage  and  allowed  direct  comparison  of 


the  data  from  different  masses. 


(3)  The  cracking,  or  fragmentation,  of  the  dimethyl  nitramine 
from  the  electron  impact  ionization,  rather  than  thermal 
decomposition,  resulted  in  substantial  increments  to  some 
of  the  lower  m/z  signals.  The  magnitude  of  this  enhance¬ 
ment  was  determined  in  two  ways. 

(a)  With  the  laser  off  (no  thermal  decomposition), 

ac  signals  (if  any)  for  each  mass  were  read  from 
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the  lock- In  amplifier.  These  readings  were  multi¬ 
plied  by  the  appropriate  gain  factor  (see  above) , 
and  divided  by  the  reading  at  m/z  90  for  the  DMNA 
molecular  ion,  to  yield  the  cracking  pattern. 

(b)  The  digitizer  was  "pre- triggered"  to  record  a  signal 
prior  to  the  laser  firing,  such  that  the  first  milli¬ 
second  or  so  of  each  temporal  profile  corresponded 
to  unheated  gas.  Subtraction  of  the  background 
resulted  in  a  value  for  the  signal  analogous  to  the 
lock-in  reading  for  each  mass.  Using  this  measure  of 
the  unheated  gas  at  each  ion  mass  a  nitramine  cracking 
pattern  was  also  obtained. 

Both  procedures  resulted  in  similar  values  for  the  cracking  pattern, 

Table  1.  The  exception  was  for  contributions  to  m/z  -44.  No  accurate  back¬ 
ground  could  be  determined  for  this  mass  since  in  order  to  remove  the  nitramine 
reservoir  the  vacuum  had  to  be  broken  briefly  with  subsequent  contamination  by 
atmospheric  CC^.  As  a  result,  only  method  (a)  could  be  used  for  mass  44. 
Improved  cell  design  is  planned  to  eliminate  this  problem. 

With  the  contribution  to  each  mass’s  signal  from  mass  90  cracking  known, 
it  was  straightforward  to  subtract  the  appropriate  fraction  of  the  mass  90 
signal  from  each  temporal  profile.  A  similar  correction  to  the  m/z  30  (N0+) 
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profile  was  required  to  take  into  account  the  cracking  of  NO2.  Separate  exper¬ 
iments  using  the  present  apparatus  to  produce  laser-heated  NO2 ,  had  given  a 
m/z  30: m/z  46  ratio  of  2.1:1.  Thus,  once  steps  1-3  above  had  been  performed 
for  m/z  46  and  30  (NO2  and  NO,  respectively),  2. lx  the  m/z  46  profile  was 
subtracted  from  the  m/z  30  profile.  The  extent  to  which  masses  other  than  90 
and  46  fragmented  could  not  be  determined  quantitatively  in  the  present  study. 


Some  contribution  at  m/z  30  is  perhaps  expected  from  the  dimethylnitroxyl 
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CRACKING  PATTERNS  FOR  DMNA  AND  NO. 


Cracking  Cracking 

(Digitizer)  (Lock-in  Amplifier) 


90 

1 

1 

74 

0.137 

0.073 

60 

0.339 

0.298 

59 

0.018 

0.027 

46 

0.126 

0.109 

45 

0.087 

0.107 

44 

-  - 

0.565 

43 

2.89 

2.02 

42 

3.36 

2.38 

30 

0.557 

0.36 

NO2  (Laser  Heated) 


radical  and/or  nitrosomethane ,  but  essentially  no  contribution  is  expected  from 
dimethylamino  radical,  dimethylamine ,  or  N-methylmethylene  imine.  In  other 
words,  it  is  highly  unlikely  that  any  first*  or  second-generation  products  of 
N-NO2  bond  split,  except  NO2  itself,  will  make  any  significant  contribution  to 
the  m/z  30  signal. 
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When  all  of  the  above  corrections  had  been  performed,  the  resulting  tem¬ 
poral  profiles  corresponded  to  changes  caused  solely  by  the  laser- induced 
thermal  decomposition  of  the  dimethyl  nitramine.  (Note  that  we  have  assumed 
the  cracking  pattern  for  dimethylnitramine  is  independent  of  whether  the  beam 
is  formed  from  laser-heated  nitramine  or  from  near  room  temperature  nitramine. 
Since  the  sampling  results  in  rotational  and  some  vibrational  cooling,  some 
"leveling"  is  expected  for  the  effect  of  temperature  on  the  cracking  pattern, 
which  are  anticipated  to  be  minor  in  the  first  place.  Furthermore ,  the  m/z  of 
most  interest  here  (30,  45,  46,  59,  and  60),  are  not  large  peaks  in  the  DMNA 
spectrum,  and  the  corrections  to  the  signal  representing  NO,  NO2,  and  the 
dimethylnitroxyl  radical  are  minor.  In  the  particular  case  of  NO2,  cracking  to 
give  m/z  30,  where  the  correction  to  the  signal  that  nominally  represents  NO  is 
large,  the  cracking  patterns  were  determined  for  laser-heated  NO  and  NO2.) 


RESULTS 


A  typical  set  of  temporal  profiles  are  shown  in  Figures  3-5.  Figure  3  is 
a  temporal  profile  of  the  buffer  gas,  argon.  One  would  expect  the  signal  level 
before  and  after  the  laser  pulse  to  remain  unchanged.  However,  the  rapid  heat¬ 
ing  of  the  reaction  mixture  propagated  a  shock  wave  within  the  reaction  zone. 

The  effects  of  this  shock  wave  are  evident  in  Figure  3.  (It  was  found  necessary 
to  machine  an  insert  for  the  aluminum  block  cell  to  fill  as  much  as  possible  of 
the  volume  unswept  by  the  laser,  in  order  to  minimize  the  effects  of  the  shock 
waves . 

These  efforts  to  minimize  the  cooling  and  reheating  that  results  from 
expansion  and  from  reflected  shock  waves  were  at  least  partially  successful, 
and  examination  of  the  m/z  40  profile  for  argon  in  Figure  3  (and  the  other 
figures,  as  well)  reveals  that  they  can  be  divided  into  several  different 
regions.  The  first  of  these  regions  is  from  zero  time  to  -  1.2  ms.  In 
Figure  3  this  corresponds  to  unheated  argon,  as  explained  above.  The  sharp 
spike  immediately  following  this  region  was  caused  by  the  electrical  discharge 
as  the  laser  fired.  It  acted  as  a  marker  for  the  arrival  of  the  laser  light  at 
the  reaction  cell.  The  large  fluctuations  in  signal  level  following  the  spike 
were  due  to  shock  waves,  caused  by  the  rapid  heating  of  the  gas,  reverberating 
around  the  cell  and  resulting  in  rapid  changes  in  pressure  and  temperature. 

When  no  SFg  was  present,  i.e.,  when  no  heating  occurred,  this  phenomenon  did 
not  occur.  The  signal  returned  to  a  steady  level  after  2.8  ms.  The  level 
during  this  quiescent  period  was  almost  exactly  the  same  as  the  level  prior  to 


the  spike.  At  -  5.0  ms  a  second  series  of  shock  waves  was  observed.  Again, 
this  phenomenon  did  not  occur  when  there  was  no  heating.  All  the  raw  data 
obtained  showed  exactly  the  same  shock  effects  at  exactly  the  same  times  with 
the  region  between  2.8  ms  and  5.0  ms  being  apparently  shock  free.  From  these 
observations,  we  concluded  that  the  region  from  2. 8 -5.0  ms  represented  the 
unperturbed  post-shock  signal  level  for  each  mass. 

The  upper  trace  in  Figure  4  is  a  temporal  profile  for  dimethyl  nitramine 
(ra/z  -  90) .  The  same  general  features  as  those  for  argon  can  clearly  be  seen. 

As  explained  above,  the  level  prior  to  the  electronic  spike  corresponds  tc 
unheated  gas  and  the  level  from  2. 8 -5.0  ms  corresponds  to  the  level  after  heating 
In  this  particular  case,  the  nitramine  underwent  -  12.5%  decomposition. 

The  temporal  profile  for  m/z  -  46,  NO2,  is  shown  in  Figure  4.  Once  again, 
the  same  general  features  can  be  observed.  The  reduction  in  intensity  of  the 
spike  and  the  different  shock  pattern  between  the  spike  and  2.8  ms  are  a  direct 
result  of  the  correction  for  parent  cracking.  The  subtractions  involved  in  the 
cracking  corrections  (see  above)  did  not  always  result  in  exact  cancellation  of 
the  spike  and  shock  patterns.  To  obtain  exact  cancellation  one  would  require 
that  the  magnitude  of  the  spike  and  shock  patterns  for  each  m/z  signal  be 
proportional  to  the  cracking  pattern  used  in  the  subtraction. 

The  shock  patterns  in  Figure  4  were  typical  for  all  m/z  signals,  once 
corrected  for  cracking.  It  must  be  reiterated  that  the  raw  data  for  each  mass 
show  the  same  shock  pattern  as  that  for  argon  and  dimethyl  nitramine.  The  zero 
level  in  Figure  4  for  m/z  -  46,  from  t  -  0  to  the  spike,  is  genuine;  that  is, 
following  background,  gain  and  cracking  corrections,  this  level  was  indeed  zero. 
This  is  as  it  should  be,  since  there  was  no  NO2  present  in  the  original  reac¬ 
tion  mixture.  Had  this  level  not  been  zero,  it  would  have  indicated  a  build-up 


of  reaction  products  between  laser  shots,  i.e.,  the  cell  contents  would  not 
have  been  fully  exchanged  between  shots.  In  fact,  if  the  gas  flow  rate  was 
decreased  or  the  laser  repetition  rate  increased,  such  a  build-up  was  observed. 
A  second  check  for  product  build-up  was  possible.  Signal  levels  for  each  m/z 
were  taken  with  the  laser  off  (no  heating)  using  a  pulse  generator  as  digitizer 
trigger  (see  above).  Signal  levels  from  t-0  to  the  spike,  obtained  with  the 
laser  firing,  were  compared  with  the  laser-off  values.  If  build-up  were  occur¬ 
ring,  one  would  expect  these  two  levels  to  differ.  No  such  difference  was 
observed.  This  was  found  to  be  the  case  for  all  other  masses,  i.e.,  the 
initial  zero  levels  are  all  genuine. 

It  is  clear  from  Figure  4  that  following  the  heating  there  is  a  sharp 
increase  in  the  NC>2  level.  In  addition,  the  rate  of  increase  of  the  NO2  cor¬ 
responds  exactly  with  the  decrease  in  the  nitramine,  indicating  that  the  NO2 
was  formed  at  the  same  time  as  the  nitramine  decomposed.  This  was  also  found 
to  be  the  case  for  the  signal  change  at  all  m/z,  with  one  exception  (see  below) 
Figure  4  thus  indicates  that  one  thermal  decomposition  pathway  for  dimethyl 
nitramine  gives  nitrogen  dioxide  as  an  initial  product. 

Since  the  nitramine  has  a  mass  of  90  and  NO2  has  a  mass  of  46,  the 
expected  counterpart  to  NO2  at  mass  46  in  the  DMNA  decomposition  is  the 
dimethyl  amino  radical,  mass  44.  There  is  only  a  very  slight  increase  in  the 
m/z  44  signal,  but  quite  clear  increases  in  the  levels  of  m/z  43  and  42, 
evidently  the  principal  (and  not  unexpected)  ion  fragments  of  mass  44. 


From  the  direct  observation  of  these  initial  decomposition  products  it  is 
clear  that  one  decomposition  pathway  involves  the  cleavage  of  the  N-N  bond  in 
dimethyl  nitramine. 

(CH3)2N-N02  -  (CH3)2N  +  N02  (1) 

However,  these  were  not  the  only  initial  products  from  the  thermal  decom¬ 
position  of  dimethyl  nitramine.  Masses  59,  45  and  30  were  also  produced 
(Figures  4-5).  Each  of  these  profiles  exhibited  the  rise  time  indistinguish¬ 
able  from  the  rise  of  m/z  46  and  the  fall  of  DMNA  at  m/z  90.  The  elapsed  time 
before  the  periphery  of  the  laser-heated  region  expands  into  the  nozzle  is  no 
more  than  4  /is;  an  upper  limit  to  the  fractional  reaction  of  dimethylamino 
radical  and  N02,  based  on  extrapolated  literature  values  for  amino  and 
methyl  radical  oxidation  by  N02,  can  be  placed  at  -  5%.  Thus,  in  order  to 
explain  the  production  of  these  species  on  a  time  scale  appropriate  only  for  a 
unimolecular  reaction,  one  must  invoke  a  nitro-nitrite  rearrangement  followed 
by  0-N  bond  scission  in  the  nitramine,  reactions  (2)  and  (3). 

(ch3)2n-no2  (CH3)2NO-NO  (2) 

<CH3)2NO-NO  -  (CH3)2NO  +  NO  (3) 

Reactions  (2)  and  (3)  explain  the  production  of  m/z  30,  NO.  Clearly,  one  would 
also  expect  to  see  the  dime thy Ini troxyl  radical  (CH3)2NO,  mass  60.  No  signifi¬ 
cant  change  was  observed  in  the  m/z  -  60  signal.  However,  the  signals  at 
m/z  59  and  45,  which  would  correspond  to  loss  of  a  hydrogen  or  a  methyl  group 
from  the  ionized  nitroxyl  species,  are  seen  to  increase  substantially  (Figure  5). 
This  fragmentation  almost  certainly  takes  place  in  the  ionizer  and,  given 
almost  no  change  in  the  mass  60  signal,  is  evidently  very  facile. 


From  Che  direct  observation  of  these  initial  reaction  products  it  is  clear 
that  a  second  pathway  for  the  thermal  decomposition  of  dimethyl  nitramine  exists 


involving  nitro-nitrite  rearrangement  followed  by  0-N  bond  scission.  Calibra¬ 
tions  using  laser-heated  NO  and  NO2  performed  under  the  same  conditions  as  the 
dimethylnitramine  decomposition,  indicate  that  the  ratio  of  rearrangement  to 
simple  bond-bond  scission  (l^/k^)  is  about  0.7. 

In  addition,  the  apparent  mass  balance  for  NO  and  NO2  production  and  DMNA 
loss  ranged  from  1.05  to  1.33.  This  was  based  on  the  NO/NO2  calibrations  and 
the  partial  pressure  of  DMNA  measured  for  the  same  set  of  flow  conditions. 
Considering  the  assumptions  made  in  determining  this  balance,  this  is  good 
agreement  and  suggests  that  the  various  spectral  subtractions,  etc.,  have  not 
introduced  any  substantial  errors. 

Only  one  other  mass  showed  any  discernable  change  after  heating.  This  was 
m/z  -  74,  corresponding  to  the  formation  of  dimethylnitrosamine ,  (CH3)2N-NO. 
Figure  6  shows  temporal  profiles  for  m/z  90,  46  and  74  obtained  under  condi¬ 
tions  where  secondary  product  buildup  was  occurring.  In  this  case,  a  very 
small  increase  in  m/z  74  occurred  following  the  laser  heating.  However,  the 
rise  time  was  much  longer  than  that  for  the  other  masses,  supporting  the  expec¬ 
tation  that  the  formation  of  dimethylnitrosamine  was  a  result  of  secondary 
reactions,  presumably  the  combination  of  (CH^^N  and  NO  radicals. 
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DISCUSSION 


The  results  presented  above  indicate  that,  in  addition  to  N-NO2  bond 
scission,  a  substantial  fraction  of  dimethylnitramine  (DMNA)  decomposes,  at 
temperatures  of  about  900  K,  by  way  of  an  intramolecular  rearrangement  to  the 
N-nitrite,  followed  by  rapid  scission  of  the  extremely  weak  NO-NO  bond.  The 
interpretation  is  straightforward  if  the  results  presented  above  are  taken  at 
face  value.  The  certainty  attached  to  these  results  rests  primarily  with  three 
factors.  The  first  of  these  is  the  extent  to  which  signal  changes  indicative 
of  chemical  reaction  can  be  separated  from  the  effects  of  reflected  shock  waves 
and  gas-den.,ity  fluctuations.  The  second  is  the  extent  to  which  spectral  sub¬ 
traction  for  m/z  30,  45,  59,  and  60  correctly  isolates  the  portions  of  these 
signals  due  to  NO  and  dimethyl  nitroxyl  radical.  The  third  factor  is  the  cer¬ 
tainty  with  which  the  observed  rise  time  of  the  30,  45,  and  59  signals  can  be 
stated  to  be  too  fast  to  be  accounted  for  by  secondary  bimolecular  oxidation  of 
dimethylamino  radical  by  NO2.  The  approach  taken  to  deal  appropriately  with 
each  of  these  factors  is  described  in  the  results  section;  these  approaches  are 
elaborated  below. 

In  considering  factor  number  (1),  fluctuations  in  intensity  that  are  due 
to  reflected  shock  waves  and  occur  after  the  arrival  of  the  first  laser-heated 
materials,  it  should  be  reiterated  that  while  there  are  short  term  (-  250  /js) 
periodic  fluctuations  in  the  m/z  40  signal  for  the  argon  carrier  due  to  density 
changes,  there  is  no  change  in  the  signal  averaged  on  a  millisecond  time  scale. 
Therefore,  changes  on  the  latter  time  scale  for  reactive  components  actually 


reflect  chemical  changes  rather  than  mere  gas-density  fluctuations.  As  stated 
in  the  results  section,  we  have  chosen  as  the  observation  window  the  interval 
when  the  gas  density  changes  reflected  in  the  m/z  40  signal  are  minimal  (3.2  to 
4.5  ms).  However,  it  is  important  to  note  that  the  results  are  not  substanti¬ 
ally  dependent  on  this  choice  of  a  time  window. 

Examination  of  Figure  4,  where  the  m/z  90,  46  and  30  intensities  are 
plotted  together,  makes  it  clear  that  there  are  three  steps  of  decreasing 
m/z  46  and  30  intensity  (2.3  to  3.0,  3.6  to  5.2  and  5.8  to  6.8  ms)  and  that 
these  are  mirrored  by  three  steps  of  increasing  intensity  in  the  m/z  90 
profile.  Thus,  decreasing  amounts  of  product  NO2  correspond  to  increasing 
amounts  of  the  starting  nitramine.  The  nozzle  is  evidently  sampling  in  these 
three  time  periods,  regions  of  successively  lower  values  of  k/dt.  While  we 
have  chosen,  for  the  reason  given  above,  to  use  the  region  from  -  3.2  to  4.5  ms 
following  the  laser  pulse,  basically  the  same  result  is  obtained  if  the  region 
just  prior  to,  or  just  after,  it  is  used.  In  these  latter  regions,  there  will 
be  somewhat  more  or  somewhat  less  nitramine  decomposition,  respectively,  and 
proportionately  more  or  less  product  appearance  from  that  decomposition.  Thus, 
while  the  fractional  decomposition  is  dependent  on  the  chosen  observation 
region,  the  distribution  of  product  masses  is  not,  making  it  very  clear  that 
the  changes  observed  as  one  samples  deeper  and  deeper  into  the  laser -heated 
region,  are  not  the  result  of  (varying  extents  of)  secondary  reaction,  but  are 
the  result  of  variations  in  temperature  and  reaction  time. 

The  second  factor  bearing  heavily  on  our  principal  conclusion  is  the  extent 
to  which  ion  fragmentation  of  species  other  than  NC>2  and  dimethylnitramine 
(which  is  quantitatively  accounted  for  in  the  spectral  subtraction)  contributes 
to  intensity  at  m/z  30,  and  thus  might  be  responsible  for  apparent  rapid 


B-19 


production  of  NO.  Dimethylamino  radical  is  not  a  likely  source  of  m/z  30, 
since  dimethylamine  itself  has  almost  no  (<  2%)  intensity  at  30.  It  is,  of 
course,  likely  that  some  of  the  m/z  30  arises  from  fragmentation  of  the  molec¬ 
ular  ion  of  che  dimethylnitroxyl  radical  (CH^^NO,  but  this  in  no  way 
invalidates  the  main  conclusion  since  the  nitroxyl  radical  Is  the  partner  of 
intramolecular ly  produced  NO.  In  fact,  it  is  not  unexpected  that  ion  fragmen¬ 
tation  of  nitroxyl  radical  (m/z  60)  might  lead  via  CH^NO  to  a  substantial 
fragment  at  m/z  30.  Accordingly,  the  intensities  in  Table  2  reflect  use  of  the 
mass  spectrum  of  nitrosomethane  itself  to  make  a  rough  estimate  of  the  contri¬ 
bution  of  the  nitroxyl  radical  to  m/z  30.  Given  that  there  are  no  other 
obvious  sources  of  mass  30,  we  are  left  with  no  alternative  but  to  assign  the 
m/z  30  intensity,  once  corrected  for  NO2 ,  (CH2)2NN02,  and  CH-jNO  cracking,  to 
unimolecular  production  of  NO. 

Finally,  confident  assignment  of  the  m/z  30  intensity  in  Figure  4  to 

intramolecular  production  of  NO  depends  on  the  ability  to  preclude  bimolecular 

reaction  of  initially  formed  dimethylamino  radical  and  NO2  as  the  source  of  the 

nitroxyl  radical  and  NO.  As  stated  in  the  results  section,  the  time  available 

for  secondary  bimolecular  reaction  is  limited  to  the  time  before  the  laser- 

heated  region  expands  the  short  distance  (<  0.1  cm)  to  the  sampling  orifice. 

With  this  expansion  taking  place  at  roughly  1/3  sonic  velocity,  the  available 

reaction  time  Is  less  than  3  ps;  using  7  x  10  lm  s  as  an  upper  limit  for 

N02  oxidation  of  dimethylamino  radical  (based  on  measured  rates  and  temperature 

dependence  for  NO2  reaction  with  methyl,  amino,  and  acetyl  radicals**'^ ,  and  the 

o 

relevent  entropy  considerations  ) ,  one  reaches  the  conclusion  that  no  more  than 
5%  of  intially  produced  NO2  could  have  been  converted  to  NO  by  bimolecular 
reaction.  Furthermore,  as  time  passes,  the  expansion  proceeds  further  and  the 
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CHANGE  IN  ION  CURRENT  (t  -  0  to  t  -  3.5  ms)  FOLLOWING  LASER  HEATING 


A  Ion  Current 


5500a 


2.9  x  10*  mol  cm* 


21310 


14100 


2.2  x  10“  mol  cm" 


Corresponds  to  ~  12.5%  depletion. 


Concentrations  determined  from  calibration  experiments  using  Ar/SF^/TM^  and 


Ar/SFg/NO  mixtures . 


Corrected  for  mass  90  cracking. 


Corrected  for  mass  90  and  mass  46  cracking. 


sampling  takes  place  deeper  into  the  laser-heated  region.  Thus,  if  secondary 


reactions  were  responsible  for  the  observed  NO,  then  the  proportion  of  NO 


should  increase  with  increasing  time.  As  described  in  the  results  section, 
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such  an  increase  is  definitely  not  observed  under  the  conditions  of  Figures  3 


through  5.  However,  if  the  DMNA  reservoir  temperature  is  raised  30°C  and  the 


laser  heating  is  increased  so  that  the  bulk  of  the  DMNA  is  decomposed,  then  the 


partial  pressures  of  initially  produced  NC>2  and  dimethylamino  radicals  increase 


by  an  order  of  magnitude  each,  and  the  rate  of  bimolecular  production  of  NO 


increases  by  two  orders  of  magnitude.  Under  these  conditions,  we  are  able  to 


see  the  NO2  signal  decrease,  the  NO  signal  grow  substantially  during  the 


observation  time  and  also  to  see  a  growth  in  m/z  74  (Figure  6) 
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CONCLUSIONS 


The  observation  of  substantial  amounts  of  NO,  together  with  the  ion  frag¬ 
ments  of  the  dimethylnitroxyl  radical,  on  the  same  time  scale  as  N0£  is  pro¬ 
duced,  and  much  too  rapidly  for  significant  secondary  bimolecular  reaction, 
necessitates  the  conclusion  that  this  rapidly  produced  NO  arises  from  an  intra¬ 
molecular  nitro-nitrite  rearrangement,  followed  by  extremely  rapid  unimolecular 
scission  of  the  very  weak  O-NO  bond.  Although  this  is  analogous  to  similar 
rearrangements  that  have  been  reported  for  nitroalkanes^  and  nitroaromatics , ^ ^ 
it  is  in  contrast  with  previous  reports  that  dime thy Ini tramine  decomposition 
proceeds  entirely  by  N-NO2  bond  scission. At  this  point,  we  presume  that 
this  reaction  proceeds,  as  those  in  the  C-NO2  cases  appear  to,  as  a  very 
"loose"  rearrangement  with  rotation  of  the  N02  not  being  substantial  until 
there  is  considerable  lengthening  of  the  N-NO2  bond.  The  nature  of  this 
rearrangement  is  now  being  explored  in  this  laboratory  (for  unsubstituted 
nitramine)  using  multiconfiguration  ab  initio  calculations. 
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Figure  1  SCHEMATIC  DRAWING  OF  THE  MACHINED  ALUMINIUM  FLOW  CELL. 

(a)  carrier  gas  inlet;  (b)  dilution  gas  inlet;  (c)  flow  cell 
pumpout;  (d)  well  for  thermocouple;  (e)  substrate  reservoir; 

(f)  retaining  ring  and  shroud  for  KC1  windows;  (g)  location  of 
nozzle;  (h)  nozzle  support;  (i)  skimmer. 

Figure  2  SCHEMATIC  DRAWING  OF  THE  MOLECULAR- BEAM- SAMPLED ,  MASS 
SPECTROMETRICALLY  DETECTED  LPHP  SYSTEM. 

For  the  sake  of  clarity,  the  diffusion  pump  associated  with  Chamber 
2  has  beem  omitted;  (a)  chopper. 

Figure  3  TEMPORAL  PROFILE  FOR  m/z  40. 

(a)  electronic  spike;  (b)  signal  fluctuations  induced  by  reflected 
shock  waves;  (c)  quiescant  period. 

Figure  4  TEMPORAL  PROFILES  FOR  m/z  90,  46,  and  30. 

Changes  in  m/z  90  are  mirrored  by  m/z  46  and  30.  The  rise  times 
for  m/z  46  and  30  are  identical  and  correspond  exactly  to  the  drop 
in  m/z  90  signal. 

Figure  5  TEMPORAL  PROFILES  FOR  m/z  45  and  59. 

Both  profiles  exhibit  the  same  rise  time  as  m/z  46  and  30. 

Figure  6  TEMPORAL  PROFILES  FOR  m/z  90,  46,  and  74. 

Obtained  under  different  conditions  than  those  in  Figures  3-5. 

The  rise  time  for  m/z  46  corresponds  to  the  fall  in  m/z  90. 

However,  the  rise  time  for  m/z  74  is  much  slower  than  for  m/z  46 
indicating  that  this  species  was  formed  by  secondary  reactions. 


Figure  2.  Schematic  drawing  of  the  molecular-beam-sampled,  mass 
spectrometrically  detected  LPHP  system. 

For  the  sake  c>  clarity,  the  diffusion  pump  associated  with 
Chamber  2  has  been  omitted;  (a)  chopper. 
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Figure  5  Temporal  profiles  for  m/z  45  and  59. 

Both  profiles  exhibit  the  same  rise  time 
as  m/z  46  and  30. 
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Figure  6.  Temporal  profiles  for  m/z  SO,  46,  and  74. 

Obtained  under  different  conditions  than  those  in 
Figures  3-5.  The  rise  time  for  m/z  46  corresponds  to 
the  fall  in  m/z  90.  However,  the  rise  tir  e  for  m/z  74 
is  much  slower  than  for  m/z  46  indicating  that  this 
species  was  formed  by  secondary  reactions. 
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ABSTRACT 


Structures  of  NH2N02,  the  four  isomers  of  NH2ONO,  dissociation  products, 
and  the  transition  state  for  nitro-nitrite  rearrangment  have  been  optimized  at 
the  MCSCF/4-31G  level.  The  wavefunction  of  the  transition  state  is 
characterized  by  two  singly-occupied  orbitals.  The  geometry  may  be  described 
as  separated  NH2  and  N02  species  held  together  by  dipole-dipole  interaction. 
Energies  have  been  obtained  by  large-scale  multireference  single  and  double 
excitation  Cl  calculations  (6-31G  basis).  The  transition  state  is  predicted 
to  lie  46.68  kcal/m  above  NH2NC>2,  0.87  kcal/m  below  the  NH2  +  N02  asymptote 
and  26.71  kcal/m  above  the  H-trans,  NONO-cis  NH2ONO  isomer,  to  which  it 
leads.  From  these  results,  isomerization  is  expected  to  be  competitive  with 
N-N02  bond  scission,  as  observed  in  recent  experiments  on  thermal 
decomposition  of  dimethylnitramine  [Stewart,  Jeffries,  Zellweger,  McMillen, 
and  Golden,  to  be  published] . 
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INTRODUCTION 


The  properties  and  decomposition  kinetics  of  compounds  containing  the  NO2 

I  O 

group  have  been  extensively  studied  both  theoretically  and  experimen- 
tally  ’  in  recent  years  because  of  the  high  energy  content  of  such 

molecules,  interest  in  unimolecular  bond  cleavage  to  form  radical  pairs  and 
the  serious  theoretical  challenges  posed  by  such  systems.  In  particular,  a 
recent  experimental  investigation^  of  the  thermal  decomposition  of 
dimethylnitramine  (DMNA)  (Cl^^NNC^  has  determined  the  rate  constant  for  DMNA 
disappearance,  log  k  (sec"'*')  -  13.5  ±  0.6  -  37.7  ±  1.5/2.3RT.  This  tempera¬ 
ture  dependence  is  not  consistent  with  the  expected  primary  decomposition 
process,  N-NO2  bond  scission,  which  has  both  a  larger  A  factor  and  a  higher 
activation  energy.  Numerical  modeling  satisfactorily  reproduces  the  observed 
temperature  dependence  and  the  observation  of  high  yields  of  dimethyl - 


nitrosamine  (DMNO)  (CH^^NNO  as  the  primary  product  only  when  rearrangement  to 
give  the  nitrite  form,  (CI^^NONO  is  included  as  a  competitive  pathway. 

Furthermore,  a  subsequent  experiment^  in  which  a  molecular  beam  sampling 
system  was  used  to  examine  directly  the  products  of  dimethylnitramine 
decomposition  found  that  the  initial  products  include  not  only  dimethylamino 
radical,  (CHj^N,  and  NO2 ,  but  also  comparable  amounts  of  NO  and  various 
fragments  ions  of  the  nitroxyl  radical,  (CH-j^NO.  The  latter  two  are  taken  as 
direct  evidence  of  nitro-nitrite  rearrangement  since  the  bond  energy  to 
release  NO  from  the  nitrite  is  expected  to  be  small  on  the  basis  of 
calculations^ ’ and  of  analogy  to  known  C-0N0  compounds. ^ 

Experimental  evidence  for  nitro-nitrite  rearrangement  in  C-nitro 

12 

compounds  is  provided  by  the  work  of  Wodtke ,  Hintsa,  and  Lee  who  applied 
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RRKM  theory  to  the  data  for  CH^NC^  to  estimate  a  barrier  height  for 
isomerization  slightly  below  the  bond  dissociation  energy. 

There  is,  to  our  knowledge,  little  theoretical  study  in  the  literature  of 

•5 

nitro-nitrite  rearrangement  of  nitramines.  Melius  has  reported  a  barrier  for 
rearrangement  of  NH2NO2  but  has  cautioned^  that  the  RHF  model  used  may  not  be 
valid  for  the  transition  state.  For  methyl  nitrate,  McKee  has  optimized  a 
transition  state  for  rearrangement  at  the  HF  level,  which  produces  a  barrier 
significantly  higher  than  the  bond  dissociation  energy.  He  conjectures  in 
this  case,  as  well,  that  use  of  a  multiconfigurational  approach  could  affect 
the  reported  results.  Although  much  of  the  published,  ab  initio  work  on  NC^- 

containing  molecules  has  relied  on  Hartree-Fock-based  methods,  it  has  been 
3  8 

pointed  out  ’  that  proper  description  even  of  the  ground  state  of  the  stable 
molecules  requires  including  more  than  one  configuration.  The  use  of  a 
multiconfigurational  approach  should  have  an  even  greater  effect  on  the 
description  of  a  transition  state. 

In  view  of  the  experimental  interest  in  decomposition  of  nitramines  and 
evidence  for  isomerization  as  well  as  the  lack  of  multiconfigurational 
treatments  for  this  important  class  of  compounds,  we  have  undertaken  ab 
initio  calculations  of  nitro-nitrite  rearrangement  of  NH2NO2 ,  the 
unsubstituted  analog  of  dimethylnitramine .  Structures  of  the  stable  isomers 
as  well  as  that  of  the  transition  state  have  been  optimized  at  the  MCSCF  level 
using  analytic  second  derivatives.  Large-scale  configuration  interaction  (Cl) 
calculations  were  performed  at  the  optimized  geometries  for  determination  of 
energy  differences.  In  the  next  section,  the  design  of  appropriate  MCSCF  and 
subsequent  Cl  calculations  is  described.  Results  and  discussion  are  presented 
in  the  following  section. 


CALCULATIONS 


Geometries  of  the  stable  isomers,  the  dissociation  products,  and  the 
transition  state  have  been  optimized  at  the  MCSCF  level  using  analytic  second 
derivatives  with  a  standard  4-31G  Gaussian  basis  set.  The  active  space  of 
the  MCSCF  calculation  was  carefully  chosen  to  include  all  the  bonds  involved 
in  the  rearrangement  of  NH2NO2  to  N^ONO.  Following  the  discussion  of 
Bauschlicher  et  al.  on  the  design  of  MCSCF  calculations  for  ^0^,  two  N-0 
bonds,  two  N-0  antibonds  and  analogously,  four  N-H  bonding  and  anti-bonding 
orbitals  were  excluded  from  the  active  space.  For  the  NH2NO2  isomer,  for 
example,  which  has  the  four  heavy  atoms  nearly  co-planar,  this  leads  to  an 
eight  orbital  active  space  consisting  of  the  four  pi  orbitals  formed  from  the 
out-of -plane  p's  on  the  heavy  atoms  and  the  following  four  sigma  orbitals:  two 
combinations  of  the  in-plane  lone  pair  p  orbitals  on  the  0's  and  the  N-N  bond 
and  anti-bond  pair.  The  eight  orbital  MCSCF  was  used  in  the  geometry 
optimization  of  the  NH2NO2  species.  However,  for  MCSCF  calculations  on  the 
NH2ONO  isomers,  a  seven  orbital  active  space  converged  far  more  readily.  In 
the  calculations  on  N^NC^,  the  lowest- lying  orbital  in  the  original  active 
space,  one  of  the  oxygen  lone  pairs,  had  an  occupation  number  extremely  close 
to  2.0.  Therefore,  removing  this  orbital  from  the  active  space  has 
essentially  no  effect  on  the  description  of  NH2NO2 .  Geometry  optimizations  of 
NH2ONO  and  the  transition  state  were  performed  with  the  seven  orbital  active 
space.  Taking  all  possible  distributions  of  10  electrons  within  seven 
orbitals  results  in  196  configurations  when  no  symmetry  is  considered. 

Energies  were  determined  by  multireference  single  and  double  excitation 
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CI  calculations  with  the  standard  6-31G  basis  set  using  MCSCF  molecular 
orbitals.  In  constructing  the  CI  expansions,  the  orbitals  were  divided  into 


three  spaces,  inactive,  composed  of  4  orbitals,  active,  composed  of  7  orbitals 
as  for  the  MCSCF,  and  external.  The  seven  lowest  occupied  orbitals  were  kept 
doubly  occupied.  For  the  stable  isomers,  which  have  Cg  symmetry,  test 
calculations  in  which  the  three  highest  frozen  orbitals  were  transferred  to 
the  inactive  space  were  performed.  For  all  Cl  calculations,  the  configuration 
list  included  all  single  and  double  excitations  from  the  inactive  and  active 
spaces  with  at  most  one  electron  in  the  external  space  and  all  double 
excitations  from  the  inactive  and  active  spaces  to  the  external  space  which 
differed  by  no  more  than  two  electrons  from  a  set  of  reference  configura¬ 
tions.  Interacting  space  restrictions  on  the  coupling  were  imposed  for  double 
excitations  from  the  inactive  space.  The  weight  of  the  reference  configura¬ 
tions  in  the  Cl  wavefunctions  was  approximately  0.88  for  all  calculations. 

This  resulted  in  Cl  expansions  of  -500,000  to  900,000  configurations  in  Cs 
symmetry  and  1.8  million  in  C^  symmetry.  Calculations  with  the  larger 
inactive  space  totalled  1  -  1.4  million  configurations  in  Cs  symmetry. 

Ol 

Preliminary  geometry  optimizations  were  carried  out  with  the  GAMESS 
program;  final  geometry  determinations  and  all  other  calculations  used  the 
ALCHEMY  11^2  program  system. 


RESULTS  AND  DISCUSSION 


Structures  of  Stable  Isomers 

Geometrical  parameters  for  NH2NO2 ,  for  the  isomers  of  N^ONO  and  for  the 
stable  dissociation  products,  optimized  at  the  MCSCF  level,  are  given  in 
Tables  I-III.  There  are  four  isomers  of  the  nitrite  depending  on  whether  the 
hydrogens  and  the  terminal  NO  are  placed  cis  or  trans  with  respect  to  the 
central  NO  bond.  The  numbering  of  the  atoms  is  illustrated  in  Figure  1.  The 
four  heavy  atoms  of  the  nitramine  are  nearly  coplanar  while  the  hydrogens  are 
bent  with  respect  to  this  plane  by  approximately  15  degrees.  For  all  four 
isomers  of  Nl^ONO,  the  four  heavy  atoms  are  strictly  coplanar  and  the 
hydrogens  are  symmetric  with  respect  to  the  plane.  There  is  a  double  bond 
between  the  terminal  nitrogen  and  oxygen  with  a  predicted  bond  length  of 
1.6 -1.7  A.  NH2ONO  is  expected  to  dissociate  readily  to  N^O  +  NO.  Comparing 
6-31G  MCSCF  energies,  the  maximum  energy  difference  between  the  isomers  is 
2.6  kcal/m.  Comparing  RHF/6-31G  geometries  and  energies,  the  isomers  are 

found  in  a  different  order  in  energy,  but  the  maximum  spread  is  similar. 

Parameters  obtained  at  the  RHF  level  with  the  6-31G  basis  as  well  as 
those  derived  from  analysis  of  microwave  spectra  for  NH2NO2  are  also  listed 
in  Table  I.  A  comparison  with  RHF  results  with  both  basis  sets  is  given  in 
Table  II  for  the  N^ONO  isomer  which  has  both  the  hydrogens  and  NONO  cis.  In 
all  cases,  qualitatively  similar  geometries  are  predicted  by  both  computa¬ 
tional  models.  Examining  the  MCSCF  and  RHF  parameters  obtained  with  the  same 
basis  set  in  Table  II,  it  is  clear  that  the  01-N2  bond  is  most  dramatically 
affected  by  the  MCSCF  model,  lengthening  from  the  RHF  prediction  of  1.416  A  to 
the  MCSCF  value  of  1.750  A.  This  is  consistent  with  the  dominant  correlation 
in  the  MCSCF  wavefunction  for  NH2ONO  being  an  excitation  to  an  01-N2 


antibonding  orbital.  It  is  possible  that  the  length  of  the  01-N2  bond  is 
overestimated  by  the  present  MCSCF  model  which  does  not  correlate  all  bonds  to 
the  same  extent.  A  modest  difference  in  geometry  due  to  basis  set  in 
calculations  at  the  RHF  level  is  also  illustrated  in  Table  II;  the  6-31G  bond 
lengths  are  all  shorter.  The  maximum  difference  is  0.07  A  for  the  01-N2  bond. 

In  the  nitramine,  the  MCSCF  model  has  the  largest  effect  on  the  N-N  bond 
length,  which  is  predicted  to  be  1.418  A  as  compared  to  1.356  A  at  the 
RHF/6-31G  level,  bringing  it  into  closer  agreement  with  the  value,  1.427  A 
deduced  from  experiment.  For  this  case,  excitation  to  the  N-N  antibonding 
orbital  is  an  important  correlation  contribution  to  the  MCSCF  wavefunction. 

In  calculations  on  NH2NO2  where  all  atoms  were  constrained  to  be  coplanar, 
MCSCF  vs  RHF  optimization  with  the  4-31G  basis  set  results  in  lengthening  of 
the  N-N  bond  by  0.07  A,  from  1.327  A  to  1.400  A.  The  total  energy  is  not  very 
sensitive  to  the  location  of  the  hydrogens.  The  MCSCF  energy  of  the  planar 
species  is  only  1.6  kcal/m  higher  than  that  of  the  optimized  geometry  in 
Table  I.  The  remaining  geometrical  parameters  are  in  reasonable  agreement 
with  those  derived  from  experiment.  Note  that  in  the  experimental  analysis, 
the  four  heavy  atoms  were  constrained  to  be  planar  and  the  N-0  bond  length  was 
assumed.  The  greatest  discrepancy  is  in  the  location  of  the  hydrogens.  The 
experimental  geometry  has  them  bent  from  the  plane  of  the  heavy  atoms  by  a 
greater  amount  and  has  a  somewhat  smaller  H-N-H  bond  angle  than  our  calculated 
values.  The  geometry  of  the  hydrogens,  however,  is  expected  to  have  little 
effect  on  the  rearrangement  process  under  investigation  here. 

Transition  State 

The  geometry  of  the  transition  state  for  rearrangement  of  NH2NO2  to  the 


NH2ONO  (H-trans,  NONO-cis)  isomer,  which  is  the  lowest  in  energy  in  these 


calculations,  is  illustrated  in  Figure  2.  The  Nl,  H2 ,  01,  and  N2  atoms  are 
nearly  coplanar  with  the  02  above  the  plane  and  the  HI  extending  below  the 
plane.  Optimized  geometrical  parameters  are  listed  in  Table  IV  and  calculated 
harmonic  frequencies  are  given  in  Table  V.  This  is  a  true  transition  state 
with  one  imaginary  frequency. 

The  geometry  may  be  described  as  separated  NH2  and  NO2  species  held 
together  by  dipole-dipole  interaction.  The  N1-N2  and  Nl-01  bonds  are 
extremely  long,  2.845  A  and  3.214  A  respectively,  and  the  remaining  relevant 
parameters  have  almost  exactly  the  values  given  for  the  isolated  species  at 
the  same  level  of  calculation  in  Table  III.  The  average  of  the  N2-01  and 
N2-02  bond  lengths  is  the  same  as  that  in  NO2  with  the  N2-01  bond  longer  by 
0.002  A  and  the  N2-02  bond  correspondingly  shorter.  As  indicated  in  Table  V, 
all  of  the  significant  frequencies  correspond  strictly  to  internal  motion  of 
one  fragment  or  the  other.  Motions  connecting  the  two  halves  of  the  molecule 
have  extremely  low  frequencies.  As  would  be  expected  for  this  loose 
combination  of  two  radical  species,  the  MCSCF  wavefunction  for  the  transition 
state  is  characterized  by  two  singly  occupied  orbitals.  This  transition  state 
could  not  have  been  located  by  a  Hartree-Fock  based  method. 

The  geometry  reported  in  Table  IV  is  identified  as  a  transition  state  by 
the  calculation  of  all  zero  gradients  and  one  negative  second  derivative.  The 
calculations  illustrated  in  Figure  3  verify  that  this  transition  state  does  in 
fact  connect  the  i'^NC^  and  N^ONO  (H-trans,  NONO-cis)  species,  as  claimed. 

The  N1-N2  bond  distance  is  taken  as  the  reaction  coordinate  for  the  left  half 
of  the  plot  describing  the  path  between  NH2NO2  and  the  transition  state. 

Fixing  the  N1-N2  bond  length  at  particular  distances  between  the  value  for 
NH2NO2  and  the  transition  state,  the  remaining  geometrical  parameters  were 
optimized  at  the  MCSCF/4-31G  level.  The  relative  energies  shown  in  Figure  3 


were  obtained  by  subsequent  MCSCF  calculations  with  the  6-31G  basis. 


Analogous  calculations  were  performed  for  the  path  between  the  transition 


state  and  the  nitrite  using  the  Nl-01  bond  length  as  the  reaction 


coordinate.  The  transition  state  is  found  to  be  the  highest  energy  point 


between  NH2NO2  and  the  H-trans,  NONO-cis  nitrite  isomer  and  to  connect 


smoothly  to  both  forms.  It  does  not,  however,  lead  to  the  H-trans,  NONO-trans 


NH2ONO  isomer.  The  barrier  to  formation  of  that  isomer  is  calculated 


(MCSCF/4-31G)  to  be  15  kcal/m  higher  than  the  transition  state.  (We  have 


considered  only  N^ONO  isomers  with  the  H's  trans  because  they  are 


approximately  trans  in  the  transition  state.) 


Energies 


Test  calculations  probing  the  effect  of  the  size  of  the  inactive  space  in 


the  Cl  expansion  were  performed  taking  advantage  of  the  C  symmetry  of  the 


stable  species.  The  energy  difference  between  one  of  the  N^ONO  isomers 
(H-trans,  NONO-trans)  and  NH2NO2  obtained  with  the  seven  orbital  inactive 


space  (denoted  CI7)  is  compared  with  that  from  the  four  orbital  inactive  space 


(CI4)  in  Table  VI.  The  Davidson  correction^  for  the  effect  of  quadruple 


excitations,  also  listed  in  Table  VI,  should  provide  a  reasonable  estimate  in 


these  calculations  because  the  weight  of  the  reference  configurations  is  very 


similar  in  all  of  the  wavefunctions .  The  energy  difference  of  21.39  kcal/m 


obtained  from  the  CI4  calculations  with  the  Davidson  correction  is  in 


reasonable  agreement  with  the  value  of  19.04  kcal/m  from  the  larger  Cl, 


justifying  the  use  of  the  four  orbital  inactive  space  for  the  transition  state 


for  which  there  is  no  symmetry  to  reduce  the  size  of  the  Cl  expansion. 


A  comparison  between  the  energy  difference  at  the  MCSCF/4-31G  and 


RHF/6-31G  optimized  geometries  is  also  provided  in  Table  VI.  At  the  same 
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level  of  calculation  (CI7) ,  the  difference  in  energy  difference  due  to 
geometry  is  only  1  kcal/m.  In  addition,  at  the  RHF  geometry,  we  compare  the 
energy  difference  from  the  present  Cl  calculations  to  that  previously  obtained 
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with  a  perturbation- theory  approach  (MP4/6-31G  ).  For  these  closed-shell 

isomers,  a  difference  in  prediction  of  3.5  kcal/m  is  found. 

The  final  results  of  our  study  are  presented  in  Table  VII  which  lists  Cl 
energies  for  the  transition  state,  the  NH2  +  NC>2  asymptote,  and  the  NH2ONO 
(H-trans,  NONO-cis)  isomer  to  which  the  transition  state  leads,  at  MCSCF/4-31G 
optimized  geometries.  The  transition  state  is  predicted  to  lie  46.68  kcal/m 
above  NH2NO2 ,  0.87  kcal/m  below  the  asymptote  and  26.71  kcal/m  above  the 
H-trans,  NONO-cis  nitrite  isomer. 

Discussion 

In  the  following  discussion,  we  assume  the  results  obtained  here  for  the 
unsubstituted  nitramine  can  be  applied  directly  to  dimethylnitramine  which  has 
been  studied  experimentally.  This  appears  to  be  a  reasonable  assumption  for 
isomerization  and  bond  dissociation  in  that  the  hydrogens  remain  as  spectators 
in  both  processes.  Because  the  frequencies  of  the  transition  state  are  very 
similar  to  those  of  the  separated  products,  consideration  of  zero  point  energy 
will  not  affect  the  conclusion  that  within  the  uncertainty  of  the  calculation, 
the  barrier  to  nitro-nitrite  rearrangement  is  the  same  as  the  energy  to  break 
the  N-NO2  bond.  Furthermore,  the  A  factors  for  isomerization  and  bond 
scission  should  also  be  quite  similar,  leading  to  the  prediction  that  the  two 
paths  should  compete  quite  evenly.  This  in  agreement  with  the  conclusion  from 
the  molecular -beam- sampled  experiment^  on  DMNA  that  the  rearrangement/N-NC^ 
bond  scission  branching  ratio  (at  ca.  900K)  is  -0.7.  While  the  present 
results  are  in  qualitative  agreement  with  the  inference  of  the  importance  of 
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nitro-nitrite  rearrangement  in  explaining  the  observed  temperature  dependence 
in  the  earlier  decomposition  study, ^  they  do  not  support  the  assumption  of  a 
small  A  factor  and  of  an  isomerization  barrier  15  kcal/m  below  the  bond 
dissociation  energy  required  by  the  kinetic  model. 

With  respect  to  the  computational  model,  our  results  indicate  the  MCSCF 
approach  is  essential  to  determining  the  open-shell  transition  state  for 
nitro-nitrite  rearrangement.  The  closed  shell  prediction  for  the  transition 
state  using  the  RHF  model'*-*’  is  characterized  by  much  shorter  bond  lengths  and 
is  significantly  higher  in  energy  than  the  transition  state  reported  here. 
Optimization  of  the  stable  isomers  at  the  MCSCF  level,  understandably,  has  the 
greatest  effect  on  the  bond  being  correlated  in  the  MCSCF.  This  effect 
appears  to  be  comparable,  for  example,  in  NH2NO2  and  in  Ni^ONO  leading  to  a 
similar  energy  difference  between  isomers  for  the  MCSCF-  and  RHF-optimized 
geometries.  Application  of  this  MCSCF  model  to  additional  N-nitro  and  C-nitro 
compounds  should  lead  to  interesting  results. 
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Table  I.  Optimized  geometrical  parameters  for  NH2NO2 
in  Angstroms  and  degrees. 


MCSCF/4 - 31Ga 

RHF/6-31G1 

N1-N2 

1.418 

1.356 

N2-01 

1.227 

1.191 

Nl-H 

0.991 

0.998 

<H-N-H 

123.3 

116.8 

C01-N2-02 

126.6 

127.0 

CH-N1-N2 

113.8 

110.5 

<Nl-N2-0 

116.7 

116.5 

<02-N2-01-Nld 

178.1 

177.6 

<Hl-Nl-N2-01d 

16.8 

25.6 

<H2-Nl-N2-01d 

164.9 

156.5 

energy  (hartrees) 
MCSCF/6-31G* 


-259.71268 


-259.71665 


Experiment 


1.42710.002 


1.206  (assumed) 


1.00510.010 


115.212. 


130.110.3 


180.0 


aEight  orbital  active  space,  including  one  additional  0  lone  pair  orbital,  was 
used  for  this  optimization. 

^C.  Melius,  Ref.  16. 

CJ.  K.  Tyler,  Ref.  23.  Angle  between  NH2  and  NNO2  plane  -  51.811. 
torsion  angles. 


(kcal/m) 

"orientation  of  H's  with  respect  to  Nl-Ol  bond. 
bC.  Melius,  Ref.  16. 
cTorsion  angles. 

^ith  respect  to  H-trans  NONO  cis. 


Table  II. 

Optimized  geometrical  parameters  of  N^ONO 
in  Angstroms  and  degrees. 

H-cisa 

NONO-cis 

isomers  < 

n 

l 

H-cis  f 

NONO-trans  5 

MCSCF/4- 31G 

RHF/4-31G 

*b 

RHF/6-31G 

MCSCF/4- 31G  ! 

1 

Nl-01 

1.384 

1.418 

1.391 

i 

1.398 

01-N2 

1.750 

1.416 

1.342 

1.720  ; 

N2-02 

1.167 

1.171 

1.159 

1.164 

Nl-H 

0.994 

0.997 

1.003 

0.995 

CH-N1-01 

113.2 

110.2 

107.0 

111.2 

<N1-01-N2 

109.9 

117.1 

117.5 

108.0 

<01-N2-02 

110.8 

114.7 

115.3 

109.6 

<H1-N1-H2 

117.0 

114.8 

115.6  $ 

<H-N1-01-N2C 

68.1 

63.8 

57.4 

65.2  “ 

<02-N2-01-N2c 

0.0 

0.0 

0.0 

180.0 

energy  (har trees) 

MCSCF/6 - 31G* 

-259.70284 

-259.70225  P 

relative  energy^ 

2.3 

2.6 

(kcal/m) 

8 

H-trans 

H-trans  K 

NONO-cis 

NONO-trans  k  & 

MCSCF/4- 3 1G 

MCSCF/4 -31G 

RHF/6-31G"  * 

Nl-01 

1.403 

1.413 

1.393  !; 

01-N2 

1.714 

1.675 

i.356  i; 

N2-02 

1.159 

1.164 

1.153  |J 

Nl-H 

0.996 

0.996 

1.002  M 

<H-N1-01 

109.5 

108.1 

104 . 5 

<N1-01-N2 

107.6 

105.5 

108.7 

<01-N2-02 

113.0 

109.1 

110.2 

<H1-N1-H2 

115.2 

114.3 

<H-N1-01-N2C 

116.4 

117.9 

123.5  1 

<02-N2-01-N2c 

0.0 

180.0 

180.0 

energy  (har trees) 

MCSCF/6 -31G* 

-259.70644 

-259.70584 

relative  energy 

0.0 

0.4 

C-17 
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Table  III.  Optimized  geometrical  parameters  for  dissociation  products 
in  Angstroms  and  degrees  (MCSCF/4- 31G) 
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Table  IV,  Optimized  geometrical  parameters  for  Nt^NOo  -*  NH2 
transition  state  in  Angstroms  and  degrees  (MCSCF/4- 31G) . 


N1-N2 

2.845 

Nl-01 

3.214 

N2-01 

1.212 

N2-02 

1.207 

Nl-H 

1.014 

<01-N2-01 

133.1 

<H2-N1-H1 

108.3 

<H2-N1-N2 

96.3 

<H1-N1-N2 

110.3 

<N1-N2-01 

96.4 

<Nl-N2-02 

111.1 

<02-N2-01-Nla 

126.9 

<Hl-Nl-N2-01a 

102.5 

<H2-Nl-N2-Ola 

9.7 

Torsion  angle. 


Table  V.  NH2NO2  -*  NH2ONO  Transition  State  Frequencies  (cm*1) 

(MCSCF/4- 31G) 


no2 

Transition  State 

nh2 

Expta 

calculation 

b  c 

Imaginary 

(184.7) 

43.7 

67.7 

91.0 

173.7 

294.1 

749.8 

752. 7d 

1319.7 

1324. 4d 

1668. 9e 

1527.6  1540 

1617.75 

1728. 0d 

3545. 4e 

3221.4  3307 

3664. le 

3311.2  3627 

lG.  Herzberg,  "Electronic  Spectra  of  Polyatomic  Molecules,"  (Van  Nostrand, 

New  York,  1966) . 

’C.  Melius,  Ref.  16. 

'R.  J.  Buenker,  M.  Peric,  S.  D.  Peyerimhoff,  and  R.  Marian,  Mol.  Phys.  43,  987 
(1981). 

lPure  NC>2  motion. 

*Pure  NH9  motion. 


Table  VI.  Energy  differences  between  NH2NO2  and  NI^ONO  (H-trans,  NONO-trans) 

6-31G  basis  except  as  noted. 


nh2ng2  nh2ono 

H-trans  NONO-trans 


total  energy 

total  energy 

(hartrees)3 

(hartrees)3 

AE  (kc 

MCSCF/4-31G  Geometry 

MCSCF 

-0.712682 

-0.705836 

4.30 

CI7b 

-1 . 283697c 

-1 . 253592d 

18.90 

CI7De 

-1.360033 

-1.329993 

19.04 

CI4b 

-1. 29359° 

-1.096175d 

20.83 

CI4De 

-1.177215 

-1.143141 

21.39 

RHF/6-31G*  Geometry 

MCSCF 

-0.716648 

-0.707247 

5.90 

CI7 

-1.284231 

-1.256609 

17.34 

CI7D 

-1.357619 

-1.328820 

18.08 

RHFf 

-0.647174 

-0.615038 

20.17 

MP4f 

-0.386522 

-0.352174 

21.56 

aWith  respect  to  -259.0  hartrees . 

bCI7  (CI4)  denotes  Cl  with  7  (4)  orbital  inactive  space  (see  text) . 

cNH2N02  Cg  symmetry  CI7:  1,080,501  configurations 

CI4:  597,551  configurations 

^NH20N0  Cg  symmetry  CI7:  1,445,944  configurations 

CI4:  816,400  configurations 

eCI7D  (CI4D)  denotes  CI7  (CI4)  with  Davidson  correction  for  quadruple 
excitations . 

^6-31G**  basis  set.  C.  Melius,  Ref.  16. 


Table  VII.  Total  energy  of  transition  state,  asymptote, 
and  NH2ONO  (H-trans,  NONO-cis) 

(6-31G*  basis  set,  MCSCF/4-31G  optimized  geometry). 


Transition  State 

nh2  + 

no2 

nh2 

(H- trans 

0N0 

NONO-cis) 

total 

energy 

total 

energy 

total 

energy 

energy3 

above  NH2NO2 

energy3  above  NH2NO2 

energy3 

above  NH2N0 

(hartrees) 

(kcal/m) 

(hartrees) 

(kcal/m) 

(hartrees) 

(kcal/m) 

MCSCF 

-0.658889 

33.77 

-0.656990 

34.96 

-0.706436 

3.92 

CI4b 

-1.050382 

49.57 

-1.049494 

50.13 

-1.098305 

19.49 

CI4Db 

-1.102852 

46.68 

-1.101460 

47.55 

-1.145406 

19.97 

aWith  respect  to  -259.0  hartrees . 
^Defined  in  Table  VI. 


FIGURE  CAPTIONS 


Figure  1 


Figure  2 

I 

Figure  3 
f 


MCSCF/4-31G  optimized  equilibrium  geometry  of  NH2NO2  and  two  NH^ONO 
isomers.  (a)  N1 ,  N2 ,  01,  02  are  nearly  co-planar.  (b)  cis  and 
trans  are  designated  with  respect  to  Nl-01  bond.  N1 ,  N2 ,  01,  02 
are  strictly  coplanar  and  bisect  ^  H1N1H2 . 

MCSCF/4-31G  optimized  transition  state  for  rearrangement  of  NH2NO2 
to  NH20N0  (H- trans,  NONO-cis) . 

Energy  (kcal/ra)  with  respect  to  NH2NO2  (MCSCF/6-31G*)  for  path 
connecting  NH2NO2 ,  transition  state,  and  Nl^ONO  (H- trans, 

NONO-cis).  Left  side:  For  fixed  N1-N2  bond  length  as  the  reaction 
coordinate,  all  other  geometrical  parameters  were  optimized.  Right 
side:  Nl-01  taken  as  reaction  coordinate. 
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Figure  1.  MCSCF/4-31G  optimized  equilibrium  geometry  of  NH2NO2  and  two  Ni^ONO 
isomers.  (a)  N1 ,  N2 ,  01,  02  are  nearly  co-planar.  (b)  cis  and 
trans  are  designated  with  respect  to  Nl-01  bond.  Nl,  N2 ,  01,  02 
are  strictly  coplanar  and  bisect  -j:  H1N1H2 . 


nh2  no2 


nh2  ono 


Transition  state 


> 

O 

C  20 


_ I _ J _ L 

15  2  30  270,2.85 

N1-N2(A)  3  20 


N1-01  (A) 


JA-8232  5 


Figure  3.  Energy  (kcal/m)  with  respect  to  NH2NO2  (HCSCF/6-31G  )  for  path 
connecting  Nl^NC^,  transition  state,  and  Nl^ONO  (H-trans, 

NONO-cis).  Left  side:  For  fixed  N1-N2  bond  length  as  the  reaction 
coordinate,  all  other  geometrical  parameters  were  optimized.  Right 
side:  Nl-01  taken  as  reaction  coordinate. 
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SYNTHESIS  OF  N , N-DIMETHYLNITRAMINE  BY 
NITRODEPHOSPHORYLATION  OF  HEXAMETHYLPHOSPHORAMIDE 

Jeffrey  C.  Bottaro,  Clifford  D.  Bedford, 
Robert  J.  Schmitt,  and  Donald  F.  McMillen 


(Submitted  for  publication  in  J.  Org.  Chem. ,  March  1988) 


SYNTHESIS  OF  N.N-DIMETHYLNITRAMINE  BY 
NITRODEPHOS PHORYLAT I ON  OF  HEXAMETHYLPHOSPHORAMIDE 

Jeffrey  C.  Bottaro,*  Clifford  D,  Bedford, 

Robert  J.  Schmitt,  and  Donald  F.  McMillen 

Energetic  Materials  Program,  SRI  International 
333  Ravenswood  Avenue,  Menlo  Park,  CA  94025,  USA. 

ABSTRACT :  Addition  of  hexamethylphosphoramide  to  anhydrous  nitric 
acid  at  0°  to  10°C  gives  a  200%  yield  of  N,N-dimethylnitramine ,  along 
with  a  12%  yield  of  N,N-dimethylnitrosamine ,  based  on  hexamethyl¬ 
phosphoramide  . 

Sir: 

1  2 

The  large-scale  synthesis  of  N ,N-dimethylnitramine  ’  poses  serious 
safety  problems  resulting  primarily  from  the  formation  of  nitrosamine 
by-product  and  the  use  of  potentially  explosive  mixtures  of  nitric  acid 

*3  / 

and  acetic  or  trifluoroacetic  anhydride.  ’  In  our  search  for  a  safe, 
pilot-plant-scale  synthesis  of  this  material  (not  currently  available 
commercially) ,  we  examined  the  ostensibly  exothermic  metathesis  of  a 
phosphoryl  amine  with  nitric  acid- -in  this  case  hexamethylphosphoramide 
(HMPA)--to  give  dime thy Ini tr amine  and  phosphoric  acid. 


This  synthesis  has  proved  viable  in  practice.  We  isolated  a  200% 


yield  of  N.N-dimethylnitramine  and  only  a  small  amount  (12%)  of  the 
carcinogenic  N,N-dimethylnitrosamine  impurity  when  the  reaction  was  run 
on  roughly  a  half -mole  scale,  based  on  starting  HMPA.  The  reaction  was 
run  between  0°  and  10*C;  ice-bath  cooling  was  required.  The  products 
were  isolated  by  neutralization  with  cooling,  followed  by  chloroform 
extraction. 

Notably,  only  two  of  the  three  dimethylamino  ligands  underwent 
nitration  to  give  a  200%  yield,  rather  than  all  three  giving  the  maximum 
possible  300%  yield  of  the  desired  product.  Further  study  will  reveal 
the  fate  of  the  currently  unaccounted-for  dimethylamino  ligand. 

The  implications  of  this  observation  for  explosive  technology  are 
significant.  A  phosphorus -based  scaffold  might  be  feasible  as  an 
incipient  polynitramine  framework,  using  variations  on  the  reaction 
described  here.  Elaborations  of  this  methodology  are  being 
investigated. 


EXPERIMENTAL 


Caution!  This  procedure  produces  traces  of  N,N-dimethylnitrosamine,  a 
known  carcinogen. 

Hexamethyl  phosphoramide  (50  g,  0.28  mole)  was  added  dropwise  to  a 
stirred  mass  of  600  g  (400  mL)  of  100%  nitric  acid  in  a  1-L  flask  cooled 
by  an  ice-bath.  The  rate  of  addition  was  carefully  regulated  to  prevent 
the  reaction  mixture  from  ever  heating  beyond  10°C;  the  addition  time 
was  approximately  90  min.  After  the  addition  was  complete,  the  reaction 


mixture  was  allowed  to  warm  to  room  temperature  over  an  additional  90 
min.  Workup  was  carried  out  by  quenching  the  lection  mixture  into  1  kg 
of  ice,  neutralizing  with  400  g  of  NaOH  with  ice-cooling  of  the  diluted 
mixture,  and  extracting  with  3  x  200  mL  of  chloroform.  Drying, 
filtering,  and  concentrating,  followed  by  crystallization  from  carbon 
tetrachloride,  gave  50  g  (200%)  of  N,N-dimethylnitramine ,  mp 
53-55°C.  The  mother  liquor  yielded  3  g  (12%)  of  dime thy Initrosamine ,  a 
yellow  liquid  with  nonequivalent  methyls  in  its  NMR  spectrum. 
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